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ABSTRACT 

We study the initial mass function (IMF) of one of the most massive Galactic star-forming regions 
NGC 3603 to answer a fundamental question in current astrophysics: is the IMF universal, or does it 
vary? Using our very deep, high angular resolution JHKsL' images obtained with NAOS-CONICA 
at the VLT at ESO, we have successfully revealed the stellar population down to the subsolar mass 
range in the core of the starburst cluster. The derived IMF of NGC 3603 is reasonably fitted by a 
single power law with index F ^ —0.74 within a mass range of 0.4 — 20 Mq, substantially flatter than 
the Salpeter-like IMF. A strong radial steepening of the IMF is observed mainly in the inner r < 30" 
field, indicating mass segregation in the cluster center. We estimate the total mass of NGC 3603 to be 
about 1.0- 1.6 X 10-^ Mq. The derived core density is > 6 x 10"* Mq\)c an order of magnitude larger 
than e.g., the Orion Nebula Cluster. The estimate of the half- mass relaxation time for solar-mass stars 
is about 10 — 40 Myr, suggesting that the intermediate- and low-mass stars have not yet been affected 
significantly by the dynamical relaxation in the cluster. The relaxation time for the high-mass stars 
can be comparable to the age of the cluster. We estimate that the stars residing outside the observed 
field cannot steepen the IMF significantly, indicating our IMF adequately describes the whole cluster. 
Analyzing thoroughly the systematic uncertainties in our IMF determination, we conclude that the 
power law index of the IMF of NGC 3603 is F = -0.741° j?. Our result thus supports the hypothesis 
of a potential top-heavy IMF in massive star-forming clusters and starbursts. 

Subject headings: stars: luminosity function, mass function — stars: pre-main sequence — stars: 
formation — HII regions — open clusters and associations: individual (NGC 3603) 



1. INTRODUCTION 

One of the most interesting properties of massive star- 
forming regions is the stellar initial mass function (IMF) . 
Since the pioneering work by Salpeter (1955), which led 
to a standard picture of the IMF, the so-called Salpeter 
IMF {dN/d\ogM cx A^^ with F = -1.35, where M is 
stellar mass), numerous efforts have been made to un- 
derstand the IMF of many types of objects such as field 
populations, stellar associations, and open and globu- 
lar clusters, covering the whole stellar mass range from 
massive OB stars down to substellar brown dwarfs. It is 
currently accepted that the IMF follows a single Salpeter- 
like power law in the high- to intermediate-mass range 
(-^ ^ 1 ^q)- It becomes flatter towards subsolar masses 
and peaks at a characteristic mass of several tenths of a 
solar mass. It then declines towards the brown dwarf 
mass range. Therefore, several analytical expressions 
have been proposed for the standard IMF, fo r example, a 
lognormal distribution (Mill er k. Scaloll979t[S calo 1986), 
a segmented power law distribution (Scalo"l998; Kroupai 
[2001) . and a combination of both (C habrier 20031 ). As 
these representations of the IMF fit reasonably well the 
various types of stellar populations, the idea of the uni- 
versality of the IMF was born. A universal IMF ba- 
sically suggests a universal star formation mechanism. 
However, it is somewhat intuitive that different physi- 
cal conditions such as the density, velocity fields, chem- 

Electronic address: yohei@mpe.mpg.de 

^ Based on observations obtained at the Very Large Telescope 
(VLT) of the European Southern Observatory (ESO) on Paranal, 
Chile, under programs 70.C-0490 and 74.C-0764 (NACO), and on 
observations from the ESO archive (ISAAC). 



ical composition, and tidal forces in the natal molecular 
clouds can lead to different star formation processes and, 
consequently, some variability in the IMF. Indeed, there 
is growing evidence for the variable IMF. There are sig- 
nificant variations in the power law index and the charac- 
teristic masses at which the distribution shows the peak 
or at which the power law breaks. A typical example of 
an IMF that does not follow the Salpeter-like distribu- 
tion is the so-called top-heavy IMF in starburst galaxies 
and young, massive star-forming regions. 

Docs the IMF really vary among stellar populations? 
To answer this question, we need to clarify, for at least 
some stellar populations, if the observed IMF variations 
are without any doubt true deviations from the Salpeter- 
like IMF or if they can simply be accommodated by the 
combined effects of observational, theoretical, and statis- 
tical uncertainties? Many studies have so far addressed 
this question. In order to address the question and to ob- 
tain new insights into the IMF of an intense starburst en- 
vironment, we study the IMF of the massive star- forming 
HII region NGC 3603 based on unprecedented spacial 
resolution observations of its central starburst cluster ob- 
tained by NACO at the Very Large Telescope (VLT), 
as well as a wider field with the Infrared Spectrometer 
and Array Camera (ISAAC) at VLT at the European 
Southern Observatory (ESO). Compared to extragalac- 
tic starbursts (e.g. M82), young star- forming regions in 
the Milky Way and in the Magellanic Clouds ~ so-called 
nearby starburst templates - are close enough to resolve 
the individual stars with current powerful adaptive optics 
(AO)-assisted ground-based telescopes and space-based 
facilities. These objects include, for example, the Arches 
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cluster, the Quintuplet cluster, R136, and NGC 3576. 
More nearby but less massive star-forming regions are the 
Trapezium cluster in the Orion Nebula Cluster (ONC), 
the Pleiades cluster, the Taurus cluster, and IC 348. 

The IMF of NGC 3603 has been presented in several 
studies, and recent works have derived somewhat flat 
IMFs with a power law i ndex of F = —0.73 for 1 — 30 
Mq in 
Mq in 
0.4- 



Eisenhau er et al.' (119981). F = -0.9 for 2.5 - 100 
-- .Sung fc B cssell (200 ¥7and F = -0.91 ± 0.15 for 
20 MfT, in fStolte et air(|2006D . As another example 
of the IMF of massive Galactic star clusters, the Arches 
cluster has shown a slightly flat IMF with F = —0.6 to 
— 1.1 for intermediate- and high-mass st ars (|Figer et al.l 
[19991: iStolte eran[2005l : iKim et al.l[200l . Studies of the 
IMF of the Arches and other clusters are presented in 

1.1. NGC 3603 

NGC 3603 is one of the most lu minous, optically visible 
HII r egions in the Milky Way jC OSS fc RadhakrishnanI 
|1969() with its global properties such as L(Hq) ~ 1.5 x 
10^9 ergs s-i (jKennicuttl [l984D, and the total mass of 
molecular clouds ~ 4 x 10^ Mq (jGrabelskv et al.lll988D . 
Its total bolometric luminosity is as large as 10 Lq, 2 or- 
ders of magnitude larger than that of the ONC, and just 
an order of magnitude smaller than the other well known 
starburst template R136 in 30 Dor in the Large Magel- 
lanic Cloud (LMC). In particular, the starburst cluster 
located in the northern part of the gigantic HII complex 
including HD 97950 - a very compact Trapezium-like 
syste m with plenty of high-mass components (i|Walbo7nl 
119731 ) - at its center, is one of the most massive and 
the densest star-forming regions known in the Galaxy. 
The starburst cluster consists of three WNL stars, six 03 
stars, and many other late O- and B-type stars in a vol- 
ume of less than a cubic light year, providingmost of the 
ionizing radiat i on in the giant HII region p>4o ffat 1983; 
ClavtonI Il986t iMoffat et all 11994 iDrissenet al., ,1995 : 



Hofmann et"aJLlll995tlCrowther fc Dessardll998[ ) It thus 
shows remarkable similarities with R136. The distance of 
NGC 3603 from the Sun has been estimated in many ear- 
lier works by means of both photometric and kinematic 
analy sis (e.g. Ivan den Berghl [l978l : IMelnick fc Grosboll 
Il982f ). and the currently accepted value is about 7 ± 
1 kpc. Owing to its intrinsic properties such as its prox- 
imity, relatively low visual extinction of only Ay =4—5 
mag, and extreme compactness and brightness, NGC 
3603 is one of the most suitable Galactic templates of 
starburst phenomena in distant galaxies. 

2. OBSERVATIONS AND DATA REDUCTIONS 

Observations were carried out in two periods, on 2003 
March 18-21 and 2005 February 14-16, using the AO sys- 
tem NACO at the VLT at ESO on Paranal, Chile. In the 
following analysis we mainly use the data obtained in the 
second observing run because of its better data quality. 
The first run data has been used for the estimate of sta- 
tistical uncertainties in our point-spread function (PSF) 
photometry. 

The central starburst cluster including the high-mass 
stellar complex HD 97950 in NGC 3603 was imaged in 
three bands with the NACO J (1.27 /xm), H (1.66 ^m), 
and Ks (2.18 /im) filters in both runs, and an additional 
band with the L' (3.80 /im) filter in the second run. The 



bright cluster core consisting of several Wolf-Rayet (WR) 
stars and OB stars was used as the AO reference source. 
The pixel scale of ~ 27 mas was used for all bands, cor- 
responding to a field of view of about 28" x 28". In the 
L'-band imaging, a very short exposure time of 0.2 s was 
chosen so as to avoid saturation due to the substantially 
high sky background level. In addition, we carried out 
short-exposure observations in JHKs bands to acquire 
photometry of the brightest members. We spent about 
one third of the whole integration time imaging a sky 
field in JHKs bands. In the L' band we spent about 
one half of the integration time. 

2.1. Data reduction 

2.1.1. NACO data 

The data reduction for the NACO data sets was per- 
formed in the standard manner of near-infrared (NIR) 
image reductions, which consists of a sky-subtraction, 
flat-fielding, bad-pixel correction, and a mosaic combina- 
tion as major procedures. We used the IRAF^ package, 
IDL, and in part the ESO eclipse package. Note that we 
have done several corrections, such as gain and pixel scale 
corrections, which were required because of the replace- 
ment of the CONICA detector between the first and the 
second data sets. We selected only frames with reason- 
able resolutions to achieve the highest resolution for the 
mosaic frames. Consequently, typical Strehl ratios mea- 
sured on individual sources around the cluster center are 
~12%, 18%, and 32% in the JHKsAi&nd mosaic frames, 
respectively. A typical Strehl ratio in the L'-band frame 
is more than 70%. 

2.1.2. ISAAC data from ttie ESO archive 

In addition to the high-resolution NACO observa- 
tions of the cluster center, we analyzed the VLT/ISAAC 
JsHKsAyaifl data sets of a wider field of view retrieved 
from the ESO archive facility. These three-band data 
have the pixel scale of 0.147" pixel^^, corresponding to a 
field of view of about 150"x 150". In addition to the basic 
reduction procedures, two particular procedures - ghost 
removal and geometrical distortion correction - were per- 
formed. Finally, a mosaic field of about 3'. 2 x 3'. 2 was 
used for the following PSF photometry. The NACO and 
ISAAC Xg-band mosaic frames are shown in Figure [TJ 

3. POINT SPREAD FUNCTION PHOTOMETRY 

We adopted the PSF photometry for the detection of 
stellar objects and the derivation of instrumental magni- 
tudes using_thB_STAPFl package implemented in 
IDL (jDiolaiti et al.ll2000[ ). This package performs a PSF 
fitting using an empirical PSF extracted from sources in 
the image itself instead of using a PSF based on analyti- 
cal models. This is particularly suitable for AO-assisted 
imaging data, which generally show a peculiar PSF, i.e., 
a sharply peaked core with widely extended wings. We 
then performed photometric calibrations using frames of 
standard stars obtained during or close to the observing 
runs. Here we only corrected a non-negligible discrep- 
ancy between the ISAAC Jg-band filter and the J-band 

^ IRAF is distributed by tlie National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation (NSF). 
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Fig. 1.— The i^s-band mosaic frames. The wide field of 3'. 2 X 3'. 2 (5.2 X 5.2 pc) with the resolution of ~ 0."46 from the ISAAC 
data (left) and the central field of 30" X 30" (0.9 X 0.9 pc) with the resolution of ~ 0".l from the NACO data (right) are shown. The 
long-exposure NACO data is scaled to enhance the faint stars; thus, the bright stars in the center are saturated. Stellar fluxes are shown 
in square-root and log scales for the NACO and the ISAAC frames, respectively. North is up and east is to the left. 



filter in the standard system. Hereafter we use J band 
for the color-corrected ISAAC Jg-band photometry. Our 
alternative calibration using 60 identical sources in the 
Two Micron All Sky Survey (2MASS)3 catalogue yielded 
fairly close zeropoints (cr < 0.04 mag). Moreover, for 
these calibrated NACO JHKs-ha,nd lists, we performed 
an additional calibration with respect to the anisopla- 
natic effect by using the ISAAC photometry lists. In 
each annulus with a step of 2" at 8" < r < 14", we com- 
puted a typical offset of mj^ACO — misAAC magnitudes 
of a source that coincides in both NACO and ISAAC and 
corrected the offsets. 

To derive detection limits in the ISAAC field of r > 
30", where the source detection is not limited by the 
crowding but simply by the detector noise, we performed 
Monte Carlo tests of artificial source detections. Defin- 
ing a 50% recovery rate as a detection limit, the limits 
are about 20.5, 19.8, and 19.7 mag in J, H, and Ks 
bands, respectively. Photometry results are summarized 
in Tab. [1] together with properties of each mosaic frame. 
In both NACO and ISAAC data the source detection is 
primarily limited by the J-band data. 

In addition, we construct a photometry list of 
JHKs{L')-ha,nd detected sources of mj < 15.5 mag in 
the innermost r < 13" region (> 95% completeness) from 
the NACO data. 

4. COMPLETENESS LIMITS 

Even in our high spatial resolution NACO data, the 
central region of HD 97950 is extremely crowded. To cor- 
rect these incompleteness effects, we h ave employed the 
empir ical method first introduced by lEisenhauer et all 
([1991 for the ADONIS observation of NGC 3603. In 
this method we use the spatial distribution and measured 
photometry of detected sources. For this we assume (1) 
the probability of the detection of a source with a given 

^ The Two Micron All Sky Survey is a joint project of the Uni- 
versity of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the NSF. 



magnitude is uniform across the whole field, (2) every 
object has a circular area around itself over which any 
sources fainter than Am mag relative to the object can- 
not be detected, and (3) at any point in the field, the 
blending effect is dominated only by the single brightest 
neighboring star. 

The technique mainly consists of three steps: (1) the 
derivation of the critical distance using the obtained pho- 
tometry, (2) the construction of the blending map us- 
ing the critical distance, and (3) the computation of 
the incompleteness correction factor based on the map. 
For this we use the J-band magnitude, which limits the 
source detections in both NACO and ISAAC data. 

1. Derivation of the critical distance 

We first computed a distribution function 
f{Amj,d). This / illustrates how often two stars 
with a difference in magnitudes Arrij and a relative 
distance d are found in the observation. Scaling 
/ by the corresponding sampling annulus, the 
density distribution _F(Ato,/, d) should converge to 
a certain value towards large radii because stars 
do not blend with each other any more. Then by 
fitting a curve to points at which the distribution 
has half the value of the plateau in each Arrij bin, 
a representative curve defining the critical distance 
can be determined. The distribution functions 
from both the NACO and ISAAC data are shown 
in Figure [21 In this analysis we also derive another 
critical distance by taking 2/3 of the value of the 
plateau in order to provide larger incompleteness 
corrections (see § I9.7P . 

2. Construction of the blending map 

Applying the derived critical distances back into 
the photometry, we then make a map. This blend- 
ing map gives for each point the limiting magni- 
tude. Thus in the map all detected sources have 
the concentric circular areas of the blending effect 
around themselves (Figure [5]). We adopted magni- 
tudes of 2MASS photometry for very bright stars 
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TABLE 1 

Summary of mosaic frame properties 



Bandl^'l 


Num. of Total integration timeW 


FWHM of PSPl'^l 


Detection 


Num. of 




framesW at the center (s) 


(arcsec) 


thresholds'^' (ct) 


detected sources 


NAOS-CONICA 


mosaics from the first run 








J (20 s X 2) 


11 440 


0.15 


7 


1807 


H (15 s X 3) 


17 765 


0.13 


7 


3163 


Ks (20 s X 2) 


16 640 


0.096 


7 


3688 


NA09-00NI0A 


mosaics from the second run 









H (15 s X 3) 
Ks (20 s X 2) 
L' (0.2 s X 150) 
Short exposure 
J (0.35 s X 30) 
H (0.35 s X 30) 
Ks (0.35 s X 30) 



127 
36 
39 
73 

15 
15 
16 



5080 
1620 
1560 
2190 

157.5 
157.5 
168 



0.15 
0.11 
0.098 
0.12 

0.16 
0.12 
0.11 



15 
10 
7 
4 

11 
10 
7 



2589 
4228 
4269 
1216 

480 
1620 
1677 



ISAAC mosaics from the ESO archive 


Js (1.772 s X 34) 36 


2160 


0.47 


17 


10106 


H (1.772 s X 34) 59 


3540 


0.49 


11 


14089 


Ks (1.772 s X 34) 58 


3480 


0.46 


8 


16105 



'^1 An integration time and a number of exposures (DIT X NDIT) for a single frame is given in the bracket. I*^' Number 
of contributing frames to the final mosaic frames after frame selections with respect to the spatial resolution.''^' Total 
integration time of the most overlapping central region. This decreases towards outer regions from the center.''*' FWHM 
of the typical PSF extracted empirically from the mosaic frame and used in the flux measurement, '"^l Detection threshold 
applied in the PSF photometry in unit of noise standard deviation derived from the mosaic frame. We chose an optimal 
threshold for each mosaic frame separately rather than using a uniform value so as to achieve an optimal detection. 



(to,/ < 10 mag) in the ISAAC data. 

3. Computation of the incompleteness correc- 
tion factors 

A completeness fraction of a magnitude bin is given 
by measuring the reciprocal proportion of the area 
over which a source with a given magnitude is de- 
tectable to the total area of the map. Then the 
inverse of the fraction is the incompleteness correc- 
tion factor. 

Completeness is not uniform but is strongly dependent 
on the stellar density across the field, and it is predomi- 
nantly associated with the radial distance from the clus- 
ter center. Therefore we applied variable correction fac- 
tors in the analysis of the radial density profile, luminos- 
ity function (LF), and IMF. The radial variation of the 
completeness is shown in Figure [31 

Finally, in order to investigate how the resulting com- 
pleteness in our method could differ from that based 
on the commonly applied Monte Carlo test of artificial 
sources, we perform such a test for the inner NACO re- 
gion of r < 13". In this test, we first generate artificial 
sources with a given magnitude using the original PSF, 
randomly distribute them onto the mosaic frame, per- 
form the PSF photometry in the same manner as the 
original photometry list, and then derive the recovery 
rate at the magnitude. We perform this procedure in 
every magnitude within mj = 12 — 21 mag. Here, we 
set the number of artificial sources at ^ 5% of the total 
number of stars in the field so as to keep low the degree of 
crowding. The resulting completeness curve is shown in 
Figure [3] ((?ray solid curve). The curve shows fairly simi- 
lar completeness to that of the inner field in our empirical 
technique (black solid curve) up to mj ~ 19 mag and a 
slight deviation at larger magnitude. One of the possible 
contributing factors to this deviation is the saturation 
of bright stars. In this recovery test, we use the NACO 



long-exposure observation, and the detection of faint ar- 
tificial sources is to some extent restricted by the bright 
outskirts of saturated stars in the field, leading to some- 
what lower completeness. In contrast, our completeness 
measurement in the empirical technique enables us to 
use both short- and long-exposure observations together, 
and hence we do not have the saturation problem in the 
measurement. Considering that we have restricted our- 
selves within the 50% completeness limit of mj ~ 19.4 
mag throughout our analysis of the LF and IMF in this 
study, we conclude that our incompleteness correction is 
in agreement with what is expected from the standard 
recovery test of artificial sources. 

5. COLOR-MAGNITUDE AND COLOR-COLOR 
DIAGRAMS 

5.1. Distance 



In the recent optical analysis ISung &: Besselll (|2004l ) 
derived a distance of 6.9 ± 0.6 kpc by fitting isochrones 
to the observed cluster sequence including the pre-main 
sequence (PMS) population in the color-m agnitude di- 
agram (CMD). Based on NIR photometrv IStolte et all 
(|2004f ) derived the distance as well as the foreground ex- 
tinction focusing on the PMS-MS transition region. They 
derived a distance modulus of (m — M)o = 13.9 mag {d = 
6.0 ± 0.3 kpc). Their photometric distance is quite simi- 
lar to the kinematic d istance of 6.1 ± 0.6 kpc reported by 
iDe Pree et al.l (|1999( ) based on radio continuum surveys 
and is slightly smaller than a nother kinematica l dista nce 
estimate of 7.7 ± 0.2 kpc bv iNiirnberger et all (|2002). 

Applying our best estimate of the isochrones (see S I5.2p 
and the interstellar extinction (§ 15. 3p . the distance mod- 
ulus of 13.9 mag provides a reasonable fit with the ob- 
served CMD. We therefore adopt {m ~ M)q = 13.9 
mag with an error of ±0.3 mag from a visual inspec- 
tion throughout the following analysis, corresponding to 
6.0 ±0.8 kpc. 
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Fig. 2. — Critical distance curve and the blending map derived from J-band photometry from the NACO (left) and ISAAC (right) data. 
Upper panels show the scaled distribution (i.e. density) F{Amj, d) in log scale and the critical distance determined based on a curve fit to 
points with half of the converging density in each Am bin. The color bar illustrates the density, and darker color corresponds to a higher 
probability to detect a star. Lower panels illustrate the blending map. The color bar indicates the limiting magnitude. 

ter constra ins the age to be not more than several miUion 
years (e.g. iMevnet fc Maeded [20M ) . Based on Huhhle 
Space Telescope (HST) observations, the three brightest 
members designated Al, B, and C in HD 97950 have 
been classified to be \ yN6ha stars, leading to the clus- 
ter a ge of < 2.5 Mvr (iMoffat et al.l [l99l lOrissen et all 
fT995HCrowther k Dessartlll998fl . Based on the isochrone 
fitting of the (Js - Ks, Js) CMD, IStolte et al.l (I200I 
suggest 1 Myr as a common age for the MS and PMS 
populations. Also, ISung fc Besselll ()2004f ) derive an age 
of 1 (±1) Myr by fitting isochrones to the optical CMD 
from the HST observations. 

To perform the isochrone fit, we create the {J— Ks, J) 
CMD of 1158 JHKs-deiecied stars in the central region 
(r < 13") of the NACO data as shown in Figure ID In 
order to fully cover the observed mass range down to 
low-mass PMS stars with masses close to the hydrogen- 
burning limit, we have used three stellar evolutionary 
models. The MS popul ation is covered by the Geneva 
stellar evolution models (jLejeune fc Schaererll200ll here- 
after LSOl). For the high-mass PMS popu lation we 
use the isochrones from iPalla fc Stahled (|1999L hereafter 
PS99). For the low-m ass PMS po pulation we use the 
isochrones from Baraffc et al.l ()199 8i hereafter B98). We 
have assumed a solar metallicity for all isochrones. 

5.2.1. Pre-Main Sequence stars 
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Fig. 3. — Completeness as a function of the observed J-band 
magnitude calculated for four distinct radial regions. The gray 
curve shows the completeness curve from the Monte Carlo analysis 
of the innermost region. 

5.2. Age and mass-luminosity relation 

In this section we describe how we estimate the age 
of the cluster and how we derive the mass-luminosity 
{A4-L) relation used later for the conversion of the ob- 
served luminosity to stellar mass. First, an upper limit 
of the cluster age can be deduced from the evolution- 
ary status of the massive stars. The presence of WR 
stars and numerous OB-type stars in the starburst clus- 
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First we compare the PMS population with the PS99 
isochrones. We restrict the comparison to the 0.3, 0.5, 1, 
and 3 Myr isochrones. Focusing on the fit in the PMS- 
MS transition region (approx. 0.7 < J — Ks < 1.4 and 
14 < J < 16 mag), we can readily choose 0.5 and 1 
Myr isochrones. However, choosing between the 0.5 and 
1 Myr is not straightforward. Therefore, we argue that 
the best- fit age is between 0.5 and 1 Myr, and in the 
following we use a 0.7 Myr isochrone constructed from a 
linear inter polation of these tw o isochrones. 

Recently IStolte et al.l (|2004D suggested a possible pres- 
ence of a sequence of equal-mass binaries, which stretches 
along the MS and PMS-MS transition regions, suggesting 
a single star formation event ~ 1 Myr ago in NGC 3603. 
We have found no obvious secondary sequence but rather 
a scattered distribution in the region. The scatter could 
be due to a combination of the photometric uncertainty 
and an intrinsic age spread of the PMS stars of several 
hundred thousand years. Alternatively, this could reflect 
a possible uncertainty in the birth line. The zero point 
in the PS99 PMS models is defined to be the birth line. 
However, the collapse and accretion processes in the pro- 
tostar formation are not fully understood, causing a sig- 
nificant uncertainty in the birth line. Thus, it is difficult 
to discuss an age spread in a young stellar population 
whose age is <1 Myr, for which the age discrepancy be- 
tween the available PMS evolutionary models becomes 
larg er partly due to dependence of the in itial conditions 
(e.g. iTout et al.l[l99l iBaraffe et aLll2002f ). 

Considering these uncertainties, we summarize that 
the age of the PMS population is around 0.5 — 1.0 Myr 
and assume an average age of 0.7 Myr. For the low- 
mass PMS stars with < 1.2 Mq, we apply the youngest 
1 Myr isochrone from the B98 model. Among currently 
available PMS evolutionary models, the B98 model has 
been reported to have the best consistency with dynam- 
ical mass determination for stellar masses below 1.2 Mq 
(|Hillenbrand fc Whitel[200l . 

5.2.2. Main sequence stars 

We compare the MS population to the 1, 1.5, 2, 2.5, 
and 3 Myr LSOl isochrones with a particular focus on 
the locations of the three bright members 1 (Al), 2 (B), 
and 3 (C) in the CMD. Figure [5] shows the CMD of the 
three stars together with the isochrones. To determine 
the age of high-mass MS stars, we need to assume ini- 
tial masses of these three stars. They are clas sified as 
WN6h-Fabs stars bv lCrowther fc DessartI (|1998l ) and the 
authors derived the current masses of about 60 — 100 Mq. 
We therefore assume initial masses of the WN stars to 
be about 80 — 120 M© (this is deduced from Figure 5 in 
the reference). Although the 2 Myr isochrone cannot be 
clearly excluded, the best match for all three stars is the 
2.5 Myr. We thus select the 2.5 Myr isochrone. 

It is worth mentioning that potential binary/multiple 
systems in the three WN stars would be insignificant 
in our age e s timat e. Based on HST NICMOS data, 
iMoffat et all ()2004[ ) suggested that Al is possibly an 
eclipsing binary system with a 30 — 90 Mq H-rich WR 
component (WN6ha) and a 25 — 50 Mq O star compan- 
ion. The sou rce B also has been suggested to be a binary 
system (e.g. iMofFatI [19831 ) . Even thoug h the source B 
is assu med to be a binary system, Crowther fc DessartI 
(|1998[ ) derived the contribution from the likely 03-5 V 



star companion of ~ 0.1 mag, being a negligible con- 
tribution in our determination. One potential uncer- 
tainty is the still unclear nature of the three WN stars. 
For example, strong hydrogen absorption features have 
been detected in their spectra, leading to the hypothe- 
sis that they are more likely hydrogen-rich, extremely 
luminous MS stars that have not evolved but have 
mimicked WR features by g enerating st rong, hot stel- 
lar winds (iDrissen et al.lll995 : Schmutz fc Drissenlll999l : 
IMoffat et al.ll2002D . IWalborn et al.l (|2002f) reported that 
they identified two 02-3.5 III(f*) stars in the core. Thus, 
if those sources are indeed such very luminous MS stars, 
the cluster's age would be to some extent smaller than 
our estimate, allowing an interpretation of a single star- 
forming event for both the PMS and MS stellar popula- 
tions. 

If our age estimate holds, we would be confronted with 
a puzzling inconsistency or a real age spread among the 
0.5-1 Myr PMS population and the 2-2.5 Myr MS 
population. One possible explanation for this is that sev- 
eral massive stars first formed in the center of the na- 
tal molecular cloud, and subsequently the strong winds 
and radiations from those massive stars have triggered 
the second generation of low-mass stars, leading to the 
star formation over several million years. Indeed ear- 
lier studies have been suggesting sequential and ongoing 
star formation from north to south throughout the gi- 
ant HII region (iMelnick et al."1989':'Hofmann et al. lll995l: 
iDe Free et "allll 999: Taoia et al. 2001; Niirnberge r et all 
I2OOI ISung fcBesselLi2004 ). Sher 25 ()Shenil96,qc iMofeU 
,1983), for example, has turned out to be an evolved post- 
red supergiant with an hourglass-sha ped nebula., whi ch is 
excited by the HD 97950 stars (Brandner et al.lll997i ). If 
these stars indeed physically belong to NGC 3603, there 
should be at least two distinct star formation episodes 
separated by ~ 10 Myr. Note that whet her or not Sher 
25 is a cluster member is not yet clear (jCrowther et al.l 
l2006h . 

In summary, since the age estimate of the MS stars is 
based merely on the isochrone fitting to the three WN 
stars and hence is still uncertain, wc conclude 2.5 Myr 
as an upper limit of the age of the cluster population. 

5.2.3. Mass-Luminosity relation 

After the selections of the best-fit isochrones for three 
stellar mass ranges, we constructed a single M-L relation 
by connecting the three isochrones so as to convert the 
observed luminosities to stellar masses for the IMF de- 
termination. Figure ini shows the constructed single M.-L 
relation in a stellar mass versus absolute magnitude plot. 
In our M-L conversion we have used all JiJiiTs-band 
magnitudes defining the stellar mass of the closest point 
on the M-L relation, for a source, in the JHKs space. 
We adopted the distance modulus of (m — M)o ~ 13.9 
mag and the uniform interstellar extinction of Ay = 4.5 
mag. 

It is worth mentioning that, although there are 
slight discrepancies between isochrones at the connect- 
ing points, these gaps are very small and thus are in- 
significant in our IMF determination. At the connecting 
point of the PMS isochrone (0.7 Myr; PS99) and the 
low-mass PMS isochrone (1 Myr; B98) at ~ 1.2 Mq, 
the gaps are (Atoj, Am^, Arnxs) ^ (0.19, 0.14, 0.14), 
fairly comparable to our photometric uncertainties at the 
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Fig. 6. — Mass-luminosity conversion relation constructed from 
the best-fit isochrones. The single relation was created combining 
the 2.5 Myr LSOl isochrone for the MS population (> 4 Mq), the 
0.7 Myr PS99 isochrone for the PMS population (~ 4 - 1.2 Mq), 
and the 1 Myr B98 isochrone for the PMS stars (< 1.2 Mq). The 
J- and Ks-hand plots are shown in the solid and dashed lines, 
respectively. 



Fig. 4.— The (J - Ks,J) CMD from 1158 J///S's-detected 
sources in the NACO field (r < 13"). The long-dashed, solid, 
and dashed curves are isochrones from the LSOl, PS99, and B98 
models, respectively. The isochrones are plotted after corrections 
for a distance modulus of {m — M)o = 13.9 mag and an average 
interstellar extinction of Ay = 4.5 mag. The arrow illustrates the 
reddening for Ay = 1 mag. The error bars on the right side give 
the photometric errors. The box indicates a zoom on the three 
brightest stars shown in Figure [5] 
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Fig. 5. — CMD of the brightest stars. The panels show the 1.5, 
2, 2.5, and 3 Myr LSOl isochrones together with the three brightest 
WR stars 1 (B), 2 (Al), and 3 (C). The designated source num- 
bers correspond to the photometry list presented in the appendix. 
Stellar initial masses are indicateci along the isochrones. 



corresponding magnitudes (e.g., error bars are shown in 
Figure |4]) . At the connecting point of the MS isochrone 
(2.5 Myr; LSOl) and PMS isochrone at - 3.7 Mq, the 
gaps are ~ 0.3 mag in all bands but still much smaller 
than a magnitude range from a corresponding mass bin 
in our IMF. These gaps are thus insignificant in our IMF 
analysis. 

5.3. Interstellar extinction 

The line-of-sight interstellar extinction towards NGC 
3603 has been estimated in previous studies, and its av- 
erage value is around E(B — V) = 1.44 mag with only 
little variation. IStolte et all (|2004f ) derived Ay = 4.0 and 
4.5 mag for the MS and PMS populations, respectively. 
Performing an isochrone fitting in the {J — Ks, J) CMD 
of the central cluster from our NACO field, we find that 
the isochrone with a shift of Ay = 4.5 mag matches best 
the MS distribution for intermediate- to high-mass stars. 
The PMS population is also well fitted with this extinc- 
tion. This can be seen in Figure[4]in which all isochrones 
are shifted by Ay — 4.5 mag. We estimate an error of 
±0.5 mag. 

In order to explore a possible radial variation of the 
interstellar extinction, we also constructed CMDs for 
different radii from the wider ISAAC field (Figure [7]). 
The stellar population in the innermost ISAAC annulus 
(7" < r < 30") reasonably resembles the central field 
from the NACO observations (Figure [4|) , and the com- 
bined isochrone fits reasonably well the entire mass range 
from the bright MS stars to faint PMS stars. Towards 
larger radii, however, we find an increasing red wards shift 
of the stellar distribution relative to the isochrone as well 
as an increasing scatter of the color distribution. This 
trend is seen up to the 55" < r < 80" annulus, and the 
distribution does not significantly change any more at 
r > 80". We find this trend in the H - Ks versus J- H 
color-color diagrams (CCDs) for the same four concentric 
annuli and conclude that the radial increase in extinction 
is A Ay ~ 2.0 mag as the maximum at the outermost re- 
gion. Previous studies have also reported a radial trend 
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Fig. 7. — Radial variations in the CMD from four radial regions 
with radii 7" < r < 30", 30" < r < 55", 55" < r < 80", and 
80" < r < 110". The solid curve shows the combined isochrone. 



of the interstellar ex t inctio n i n NGC 3603 ( Mclnic k et al.l 
[l98l IPandev et all [200l . ISung fc BesselL (.2001 de- 
rived E{B — V) = 1.25 mag at the core (r < 0'.2) and an 
increase of 1.8 mag or even higher at larger radii. The 
radial trend can be explained by the strong stellar winds 
and radiation from the central high-mass stars which cre- 
ate a cavity in the inter stellar matter, as has been sug- 
gested in earlier stu dies (jFrogel et al.lll977l ; iBalick et al.1 
119801: IClavtol]|[l986h . 

5.4. Infrared excess emission and disks 

We construct a Ks — L' versus J — H CCD for the cen- 
tral region to examine the fraction of stars with infrared- 
excess emissions. The L-band observation is a par- 
ticularly effective tool t o identify the excess emission 
from circumstellar disks (|Lada fc Adams! [l992f ). The L- 
band observation has been used for tracing stars with 
circumstellar disks in star-formin g regions such as the 
Taur us- Auriga molecular cloud (iKenvon fc HartmannI 
|1995[ ). the Trapezium cluste r (iLada et al.l l2000l). th e 
Chamaele on I dark cloud (jKenvon fc Gomed |200TI ) , 
30 Dor (iM aerck erfc Burtmll 120051). and NGC 3576 
(|Maercker et al. 20"06|). 

Figure [8] shows the CCD for 369 stars detected simul- 
taneously in the JHKsL' bands in the cluster center 
of r < 13" after the rejection of field stars (see § 16. ip . 
Because of our shallow L'-band detection, we set the 
L'-band limit of 15.2 mag below which the detection is 
reliable. Thus, this analysis is restricted to stars more 
massive than ~ 1.2 Mq. To classify stars with and with- 
out the i'-band excess emission, we define a threshold 
of the Ks — L' color by adopting the reddening law and 
a typical error bar of Ks — L' = 0.20 mag, which is the 
1 cr scatter of Ks — L' color for the bulk of the PMS 
population. 

In total 89 out of 369 sources are classified as stars 
with disk emission, corresponding to a disk fraction of 
~ 24%. The resulting disk fraction, however, strongly 
depends on the adopted selection criterion. If we adopt 
a zero threshold instead, the disk fraction rises to ~ 38%. 
If we use a 2 (7 threshold, then the outcome is ^ 13%. 



Defining these values as a rough error bar, we conclude 
that the disk fraction for the given population and field 
is approximately 24 ± 10%. 

Parenthetically, we note that, since there are potential 
star-disk systems that are detectable only at longer wave- 
lengths, e.g., at the KsL' bands or merely the L' band, 
the use of stars detected in JHKsL' bands simultane- 
ously could underestimate the actual disk fraction. To 
roughly test this uncertainty, we compute a disk fraction 
using also stars that are detected only in L' band simply 
assuming that such stars are star-disk systems. Com- 
bining the L'-band-only detected stars and the JHKsL' 
stars with excess emissions increases the disk fraction by 
approximately 20%. 

In addition, we examine a mass dependency of the 
disk fraction. The disk fraction in stellar clusters have 
shown some dependencies on the stellar mass or spec- 
tral types. Based on N-body dynamical simulations and 
simul ations of the mass lo ss from the star-disk encoun- 
ters, iPfalzner et al.l ()200 6') find that in the Trapezium 
cluster the circumstellar disks around high-mass stars 
dissipate much faster than those of intermediate-mass 
stars, largely due to gravitational interactions. This 
shorter disk dissipation timescale for higher-mass stars 
is generally consistent with most observations. Based 
on the JHKL o bservations of the Trapezium cluster, 
iLada et "aTl (|2000f ) report that they find a considerably 
lower disk fraction for the high-mass O, B, and A stars 
than that for F-M type stars, and this could be due to 
either a lower probability for disk formation or more 
rapid disk dispersal times. Recently, based on NIR 
and mid-IR data taken by the Spitzer Space Telescope, 
iHernandez et all ()2007D study the 3 Myr old a Orionis 
cluster, and derive a disk fraction of ^ 30% for T Tauri 
stars with < 1 Mq and a much lower disk fraction for 
higher mass stars (e.g., ~ 15± 7% for Herbig AeBe stars 
with > 2 Mq), implying a more rapid disk evolution for 
stars with > 1 Mq in the cluster. 

Considering the low-mass limit of ^ 1.2 Mq imposed 
by the i'-band detection limit, we here compute the disk 
fraction for high-mass MS stars and intermediate-mass 
PMS stars so as to test the mass dependency of disk 
fraction in this cluster. Setting the MS turnon mass of 
~ 4 Mq as the threshold, we derived the disk fraction of 
^ 22% for stars with > 4 Mgand ~ 27% for stars with 
~ 1.2 — 4 Mq. As is the case of the disk fraction of whole 
mass range, we derived an error bar of about ±10% for 
the both fractions. 

Our result is consistent with earlier studies of NGC 
3603, as well as with the general consensus on the circum- 
stellar disk frequencies in intense star-forming regions. 
Based on th e L^-b and luminosity and the Ha emission, 
IStolte et all (|2004D derived a disk fraction of 20 ± 3% for 
all stars in the inner region (r < 20") and that of 12% for 
MS stars (0-A stars) only. They also reported a radial 
increase of disk fraction to ~ 40% in the outer regions 
(r < 33"). As our data is restricted within the innermost 
r < 13" region, our result is still compatible with the hy- 
pothesis that the disk fraction increases to some extent 
towards outer regions in the cluster. 

Based on a s tudy of evolutionary m odels for proto- 
planetary disks, lArmitage et al.l (2003) reported that in 
young clusters up to 30% of stars lose their disks within 
1 Myr, the remainder, however, have disk lifetimes of 
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Fig. 8. — The Ks — L' vs. J — H diagram from stars si- 
multaneously detected in the JHKsL' bands in the cluster core 
(r < 13"). The colors are dereddened using the extinction law 
bv lRieke fc Lebofs kv (1985 ). The solid curve shows the empirical 
colors of MS stars from B cssell &: BrettI II1988I) with spectral types 
from B8 to MO. The dashed line indicates the reddening vector 
passing through the B8 star. The dotted line indicates the selection 
criterion for the disk population, which is the reddening lino of the 
B8 star shifted by the 1 cr photometric error of A(Kc; — L') = 0.20 
mag. Crosses are stars with disk emissions. 



typical ly in the range o f 1 — 1 Myr. In their study of 
cr Ori. iHernandez et al.l (|2007l . Fig. 14) compile the disk 
fractions of T Tauri stars, as a function of age, for dif- 
ferent stellar populations, showing a clear correlation be- 
tween the disk fraction and the mean age. Considering 
the age of NGC 3603 of < 2.5 Myr, our disk fraction of 
approx. 30% for intermediate-mass PMS stars (1.2 — 4 
Mq) appears be somewhat smaller than those of other 
clusters, 40 — 60% based on the correlation. One reason 
could be that our analysis is limited for PMS stars with 
masses down to ~ 1.2 Mq, and the disk fraction is ex- 
pected to be much higher for lower mass T Tauri stars. 
Another reason could be that our analysis i s restricted 
to the central r < 13". As mentioned above, Stolt e et al.l 
(|2004*) report a radial increase of the disk fraction up to 
~ 40%. 

6. LUMINOSITY FUNCTIONS OF NGC 3603 

The LF described as ip{M) = dN/dM is one of the 
most important diagnostic tools in the study of star clus- 
ters, since it is not influenced by any theory-related un- 
certainties, but is purely an observational quantity. In 
contrast, the IMF relies on the transformation of the ob- 
served LF using a M-L relation and in that sense imposes 
additional assumptions. The relation between the IMF 
£,{A4) and the LF of a stellar population with a given age 
is 



ip{M) cx C(logA^) 



1 dM 

mIm' 



(1) 



Since the stellar population in NGC 3603 is young 
enough not to be strongly affected by stellar evolution, 
it provides an empirical basis for the derivation of the 
initial mass distribution. Here we note that, since we 
construct the LFs of NGC 3603 by combining the NACO 
and ISAAC observations, we verify that this combination 
is feasible by comparing two LFs for an overlap region 
at r = 10" - 13" in both the NACO and ISAAC fields. 




0.0001 Ll 



Fig. 9. — Comparison of /C^-band LFs based on three different 
estimates for the field star subtraction. The observed Kgl^F for 
9158 stars detected simultaneously in the JHKg bands within r < 
110" is subtracted by the field LFs derived from the color cut [black 
line), the outermost region (80" < r < 120"; dashed line), and the 
NPG02 KhF {gray solid line). 

We also checked that the two data sets are compatible 
for the following IMF construction. 

6.1. Cluster membership 

To derive the LF of the cluster stars, we need to cor- 
rect the potential contribution from field stars. For that 
we examine two different technical approaches. The first 
approach is to directly measure the density and LF of 
the field stars. This is a commonly applied method in 
the analysis of stellar clusters and requires the observa- 
tion of a so-called control field. Since we do not have 
a control field in our observations, we derive an upper 
limit for the field star density from the outermost region 
(80" < r < 120") of the ISAAC field. Another estimate 
is available from observations of NG C 3603 presented 
in lNiirnberger fc Petr-GotzensI ()2002l hereafter NPG02). 
Their field KLF is derived from a field out of the clus- 
ter and is in agreement with estimates based on Galactic 
disk models. The second approach is to directly differ- 
entiate the likely non-cluster members from the cluster 
members on the basis of their locations in the CMD. We 
hence call this method color cut hereafter. 

Figure [9] compares the i^s-band LF after the correc- 
tion based on the three estimates - from the outermost 
region, the NPG02 field KLF, and from the color cut. 
We first construct an observed KhF for stars simultane- 
ously detected in the JHKs bands in the whole field and 
then subtract the K\jF by the three field KhFs. We first 
find that the outermost field i^TLF overestimates the field 
star density. This is most likely due to the fact that the 
outermost region in our ISAAC field is not far enough 
and is still within the cluster region. While the T^TLFs 
subtracted by the color cut KFF and the NPG02 field 
KhF show fairly similar distributions within the reliable 
magnitude range of rriKs 12—18 mag. This verifies 
that, despite of the lack of control fields, the color cut 
method is a reasonable treatment for the field star cor- 
rection. We thus use the color cut in the following LF 
and IMF analysis. 

The selection of the cluster stars by the color cut is 
illustrated in the CMD in Figure [TOl The color cut {gray 
line) is created by smoothing out the turnover feature in 
the PMS-MS transition of the combined isochrone and 
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shifting this curve blueward by A( J — Ks) ~ 0.26 mag. 
Note that the shift is applied only in the J — Ks direc- 
tion, i.e. no shift in the vertical J direction. In total 
7514 out of 9158 stars are found to be cluster members. 
Likely noncluster stars that are located around the MS 
(dashed curve) at nij > 16 mag and bluer than the color 
cut could be part of the following groups: (1) background 
MS stars with Ay — 4.5 mag, (2) foreground low-mass 
MS stars, and (3) cluster stars formed in earlier gener- 
ations. Thus, they are excluded from the LF and IMF 
derivations. Sources in the first group would be MS stars 
with an interstellar extinction similar to those of the clus- 
ter members, but with a smaller brightness due to their 
larger distance. In this case the interstellar material such 
as gas and dust are assumed to exist locally along the line 
of sight, and hence the interstellar extinction behaves as a 
step function rather than a linear increase with distance. 
Some supposedly noncluster stars could indeed be cluster 
members that have formed prior to the starburst ^2.5 
Myr ago when the majority of the cluster members were 
born. Since the IMF study is basically the subject of a 
single star formation, such sources are correctly rejected 
from our analysis. In fact, the number of stars from ear- 
lier generations is thought to be very small considering 
the lack of sources in the low-mass MS regime in optical 
observations. Sung & Bessell (2004) found no MS popu- 
lation for M < 4 Mq in their optical observatio n. The 
absence of low-mass MS stars is also reported in iGrebell 
(|2004). 

One potential limitation of this color cut for the field 
star estimate is that the color cut curve is placed based 
on an eye inspection; therefore, the correction essentially 
relies on this judgment. However, as shown in Figure fTOl 
(and also Figure [7|) , there is a somewhat visible gap be- 
tween the low-mass PMS population and the faint MS 
population, and the color cut curve is satisfactory in sep- 
arating them. Another limitation is the fact that the 
region redward of the color cut could still contain some 
intrinsically very red non-cluster objects such as the line- 
of-sight red giants. However, the above analysis of the 
NPG02 field KLF and our KLF based on the color cut 
suggests that the contamination of red giants is statisti- 
cally insignificant in our study. Nevertheless, we study 
the potential effect on the IMF applying a "red cut" in 
addition to our color cut in § 19.81 

6.2. The cluster star LF 

Figure El shows the KLF of NGC 3603 applying the 
color cut and the incompleteness correction. The result- 
ing LF shows a monotonic increase towards faint magni- 
tude with a power law index of a = 0.27 ± 0.01. Because 
of the strong variation of the detection limit across the 
field, we cannot conclude if the observed sign of flatten- 
ing of the LF at m^s ~ 17.5— 18 mag is a real feature or 
a result of the observational limitation. We find that the 
slope of our LF is consistent w ith those of the pre vious 
studies presented in NPG02 and' Brandl et al.l ()19990 . but 
our high angular resolution NACO data reveals ~ 50% 
more stars in the central region of r < 33". 

7. INITIAL MASS FUNGTIONS OF NGC 3603 

7.1. Initial Mass Function determination 

Applying the color cut and the incompleteness correc- 
tion, we derive the IMF for the 7514 JHKs-detected 
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Fig. 10.— The (J - Ks, J) CMD for 9158 stars in the whole 
field r < 110". The black solid curve is the combined isochrone. 
The gray curve illustrates the color cut created by smoothing 
out the turnover feature in the PMS-MS transition of the com- 
bined isochrone. The dashed curves represent the extended MS 
isochrones. The non-cluster members based on the color cut are 
shown in gray squares. 
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Fig. 11.— The Xs-band LF of NGC 3603 derived from 7514 
J H K g-detected stars in the entire field (r < 110") after the color 
cut. The circles and squares show the distribution with and with- 
out applying the incompleteness correction, respectively. The verti- 
cal dashed line represents the 50% completeness limit of ~ 17.4 
mag (derived from mj ~ 19.4 mag and typically J — Kg = 2 mag 
for low-mass PMS stars) within r < 30". The vertical dotted line 
indicates the detection limit in the outer ISAAC field derived from 
the artificial source detection test. The best-fit power law slope 
{solid line) is measured in a magnitude range from 12 mag up to 
the 50% completeness limit. 

stars in the whole field. The resulting IMF is shown in 
Figure [121 To compute the stellar masses, we used the 
combined A4-L relation created from the three best-fit 
isochrones. As mentioned in § 1, the IMF is generally 
known to follow a power law ^(logA^) cx A^'", or alter- 
natively expressed in the number of stars per unit mass 
interval. 
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Fig. 12.— IMF of NGC 3603 derived from 7514 cluster stars in 
the whole field extending out to r 110". A double size bin is used 
at around 4 Mq to smooth out a discontinuity from connecting the 
MS and PMS isochrones. Diamonds and circles indicate the raw 
and the incompleteness corrected mass distributions, respectively. 
The best-fit power law slopes, derived within 0.4 — 20 Mq, are 
shown as a dotted and a dashed line. The Poisson counting error 
is used in a weighted linear least-square fit. The vertical dashed and 
dotted lines indicate the 50% completeness limit within r ~ 30" 
and the detection limit in the ISAAC outer fields, respectively. 



thus, dN/dM cx , with 7 = T - 1. The slope is de- 
rived for the mass range of 0.4—20 Mq. This mass range 
is not affected by the saturation in the ISAAC data, and 
it guarantees a completeness of at least 50% even in the 
most crowded central region. The best-fit power law in- 
dex is r = —0.74 (±0.02). The error is merely a formal 
fit error. The real errors of the index involves many other 
systematic uncertainties, and we discuss them in detail in 
§ini Although we find a turnover at around 0.2 Mq, we 
can not answer if the turnover is intrinsic or it is simply 
caused by the detection limit. 

In summary, our resulting IMF of NGC 3603 follows 
a power law with F -0.74 within 0.4 - 20 Mq. The 
detailed discussion in §[S]will show if this finding is signif- 
icant within the systematic uncertainties. In particular 
we will test if NGC 3603 as a whole or merely the ob- 
served region within r < 110" has a shallow IMF. 

7.2. Radial variation of the IMF 

One of the characteristics of NGC 3603 that makes 
it a particularly interesting object is the high concen- 
tration of massive stars in the central starburst cluster. 
This brings up the question of a possibly mass-segregated 
stellar distribution in the cluster. Utilizing the richness 
of our data, to answer the question, we have built IMFs 
of seven concentric annuli with radii of 5", 10", and 13" 
(from NACO) and 30", 55", 80", and 110" (from ISAAC) 
as shown in Figure [T51 

Here we note that there is some indication of a broken 
power law at high masses for the inner three regions from 
the NACO data in which there is no saturation problem. 
We, however, adopt the single power law with the fixed 
mass range of 4 — 20 Mq to keep consistency and to 
trace the radial variation. The IMF power law indices are 
summarized in Tab. [5) There is an obvious steepening 
of IMF with increasing radius up to r '--^ 30". For larger 
radii the slope stays more or less constant around F ~ 
—0.8. This is a clear indication of mass segregation in the 
cluster. The main characteristics of this mass segregation 
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Fig. 13. — IMFs for seven concentric annuli with increasing 
radius. Th e symbols and lines are identical to the global IMF in 
Figure [T2I 



TABLE 2 

Radial variation of the IMF power law index 



Regions 


r 


NACO data 


r < 5" 


-0.31 


5" < r < 10" 


-0.55 


10" < r < 13" 


-0.72 


ISAAC data 


13" < r < 30" 


-0.75 


30" < r < 55" 


-0.80 


55" < r < 80" 


-0.86 


80" <r < 110" 


-0.83 



are as follows: (1) a strong concentration of high-mass 
stars in the very center at r < 13", (2) a shallow IMF 
(F ~ —0.8) at r > 30", and (3) no evidence for further 
steepening at larger radii. The third point allows us to 
conclude that the IMF of the whole cluster (including 
the regions not covered by our observations) cannot be 
steeper than an IMF with F ^ -0.9. 

7.3. Comparison with previous studies 

The IMF of NGC 3603 has been studied in earlier 
works based on both ground- and space-based observa- 
tions. Based on ggT/Planetary Camera 1 (PCI) obser- 
vations jMo^IeLalJ (|1994l) derived the IMF of high-mass 
stars with a power law index of F = — 1.4zfc 0.6 wi thin 
a mass range of 30 - 60 Mq. iHofmann et all (fl995h de- 
rived the IMF with F 1.59 for 15 - 50 Mq based on 

their diffraction-limited speckle observations. These re- 
sults are consistent with the standard Salpeter IMF; how- 
ever, only the high-mass population is covered in these 
studies. 
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Here we compare our results with several recent works. 
Based on high-resolution NIR observations using the AO 
system ADONIS/SHARPII+ at the ESO 3.6 m tele- 
scope, lEisenhauer et al.l ()1998f ) were the first to reveal 
the PMS population with masses as low as ^ 1 Mq in 
the central starburst cluster. They derive a shallow IMF 
with r ~ —0.73 within 1 — 30 AfQand find no obvious 
turnover in the IMF. Our resulting IMF with F = -0.74 
over 0.4 — 20 Mq is fairly similar to their result. Even 
though we can probe much fainter stars, we still see no 
sign of turnover. 

B ased on the seeing-limited UBVRI CCD photome- 
try, [SagM^eFian (|2001 ) derived the MF for seven distant 
open star clusters, including NGC 3603. The slope of the 

MF is F 0.85 for 7-75 Mq. Ahhough then IMF is 

limited for high-mass population, our result is in agree- 
ment with their result. The authors, however, regard it 
as being consistent with the Salpeter IMF within errors. 
They mentioned that, above 1 Mq, MF slopes of star 
clusters younger than 500 Myr (equivalent to the IMF) 
in the solar neighborhood had no dependence on Galactic 
longitude, Galactocentric distance, and cluster age, and 
they were in agreement with the Salpete r IMF within 
a 1 CT error (err ~ 0.3). As mentioned in IScalol (|2005D . 
this implies that their actual uncertainties are very large. 
In contrast, our interpretation of the resulting slope is a 
case of a flat IMF. 

Recently, combining multi- wavelengt h ground-based, 
HST and Chandra X-ray observations, ISung &: Besselll 
(j2004[ ) derive a moderately flat IMF with F = -0.9 for 
2.5 — 100 Mq. The IMF shows a gradual steepening to- 
wards the outer regions (F — —0.5 at r < O'.l, —0.8 at 
r = O'.l ~ 0'.2, -1.2 at r > 0'.2). Like us they find 
a radial variation in the IMF. Our IMFs, however, are 
somewhat flatter than their IMFs. Their IMF of the out- 
ermost region (r > 12") shows a slope of F = —1.2, which 
is ~ 0.3 steeper than that of the whole field (F — —0.9), 
likely due to the mass segregation. In contrast, our IMF 
of the outer region (F ~ —0.85) is not substantially 
steeper than the IMF of the whole cluster (F ~ —0.74). 
As for the innermost region, our IMF with F ^ —0.3 for 
r < 5" is slightly flatter than then IMF with F -0.5 
for r < 6" . Th i s mig ht be in part due to the fact that 
ISung fc Besselll (|2004[ ) includes stars with masses up to 
100 Mq, while our analysis is restricted to masses below 
20 Mq. Moreover, our mass range covers down to 0.4 
Mq, which is slightly lower than their low- mass cover- 
age. 

Most recently, based on the JHKsL' ISAAC data and 
Hq, imaging fro m the HST /Wide F ield Planetary Cam- 
era 2 fWFPC2l. [Stolte et all (l200l derived the MF tak- 
ing into account the held star population, individual red- 
dening, and the potential binary contribution in the cen- 
tral region (7" < r < 65"). The derived IMF has a 

power law slope of F 0.91 ± 0.15 for 0.4 - 20 Mq. 

Compared to, for example, their result of F = —0.87 
for 7" < r < 33", our IMF shows a slightly shallower 
power law F = —0.75 for 13" < r < 30". Since the same 
data set is used, this difference could arise from techni- 
cal differences such as the field star subtraction and the 
incompleteness correction. Indeed, it appears that we 
applied a slightly smaller incompleteness correction than 
theirs. Although we present it in detail below (ii l9.7|) . we 
note that an application of a stronger incompleteness cor- 



rection is found to steepen the IMF with AF ~ —0.13. 
Thus, some part of the slight difference between both 
studies would be ex plained by the diffe rence of the cor- 
rection rates. While lStolte et all ()2006D claim no signifi- 
cant variation in power law index of the IMF but talk of 
a depletion of the high-mass tail of the stellar mass dis- 
tribution with increasing radial distance, we see a strong 
variation with radial distance from F ^ —0.3 to ~ —0.8 
over the mass range of 0.4 — 20 Mq. A part of the reason 
for this discrepancy is that Stolt e et al.l ((2006') could not 
probe the central r < 7" because of crowding, while our 
high resolution NACO observations resolve the central 
cluster even at r < 5". 

In summary, we find that our results are partly in 
agreement with previous studies, but the IMF shows a 
fairly shallower slope with F = —0.74 when compared 
to the Salpeter-like IMF. However, to prove conclusively 
the shallowness of the IMF, we have to first analyze the 
observed mass segregation. 

8. MASS SEGREGATION IN NGC 3603 

The mass segregation in a stellar cluster can potentially 
lead to a systematic error in the IMF determination if the 
observation does not cover the entire cluster. This is be- 
cause higher-mass stars are more concentrated towards 
the cluster center than lower mass stars. Since we find 
the mass segregation in NGC 3603 based on the radial 
steepening of the IMF, we need to investigate whether 
the mass segregation is the cause of the shallow IMF. 
Therefore, in this section, we study the characteristics of 
the observed mass segregation in NGC 3603. For that 
we investigate the dynamical evolutionary state of the 
cluster by estimating the relaxation time of the cluster 
(eq. [8] in 18. 3p . We first derive the global properties of 
the cluster such as its size, total mass, core radius, and 
the half-mass radius. These quantities allow us to esti- 
mate the relaxation time and then to assess the potential 
impact of the mass segregation on the IMF determina- 
tion. 

8.1. Radial mass density profile 

To determine the global parameters, we measure the 
projected radial mass density distribution and fit analyti- 
cal models to the profile. First we fit the pro file using the 
standard empirical formula bv iKinej (|1962f l. This King 
model is commonly used to describe the radial stellar 
density and brightness density profiles of globular clus- 
ters and old Galactic clusters. It is described as 



f{r) = k 



1 



1 



(3) 



where fc is a normalization factor, which is approximately 
the central mass density; rc is the core radius, and rt is 
the tidal radius. Because our color cut may not exclude 
all field stars, we modify the King model by adding a 
constant term. 

For comparison, we also fit a modified power law used 
by Elson, FaU, & Freeman (1987; hereafter EFF87) in 
their study of the radial surface brightness profiles of 
young clusters in the LMC. They find that the young 
LMC clusters do not appear to be tidally truncated even 
at radii of several hundred arcsec. The EFF87 model is 



The IMF of NGC 3603 



13 



10.00 



1.00 



0.10 



0.01 





k ~ 4.2 




- 4". 8 






N 


V BG ~ 0.056 









1 10 100 

Distance from the cluster center [orcsec] 

Fig. 14. — Projected radial mass density profile of NGC 3603 for 
intermediate- and low-mass stars within the mass range of 0.5 — 2.5 
Mq. The mass density of the innermost (r < 13") and outer field 
(13" < r < 110") is measured using 1" and 8" steps, respectively. 
The color cut and the incompleteness correction are applied. The 
solid and dashed curves show the best-fit King and EFF models, 
respectively. The horizontal dotted line shows the background level 
from the King model. 



described as 



fir) = /o 1 + - 



-7/2 



(4) 



where /o is the central surface mass density, a is a mea- 
sure of the core radius, and 7 is the power law index at 
large radii. The parameter a is related to the core radius 

of the equivalent King model by Vc = ai^l^ — 1)^/^. 

Figure [13] shows the radial mass density profile of NGC 
3603 for stars within the mass range of 0.5 — 2.5 Af©, 
which corresponds to about m^s = 15 — 17 mag. Note 
that we use the barycenter, which we derive from the 
stellar density distribution, as the center of the cluster. 
In the fit we omitted the innermost r < 2" where the 
crowding effect cannot be corrected any more. We select 
this mass range considering the following aspects. In a 
mass-segregated stellar cluster, the shape of the mass 
density profile and consequently its characteristic radii 
depend on the considered stellar mass range. Since our 
eventual goal in this analysis is to determine how many 
intermediate- and low-mass stars potentially reside out- 
side the observed field of view, we chose the upper mass 
limit of 2.5 Mq. The lower mass limit of 0.5 Mq is set 
to keep > 50% completeness. From the fit of the King 
model we derive a core radius of ~ 4". 8 (~ 0.14 pc 
at d ~ 6 kpc). Since the power law index in the best- 
fit EFF model is 7 = 1.97, it is essentially iden tical to 
the King model with — > 00. For comparison, iGrebell 
([2004) report a core radius of Tc ~ 0.25 pc based on RST 
observ ations, and NPG02 re port ^ 23" as an upper 
hmit. ISung fc BesselH ()2004f ) derive - 3" from their 
optical data. 

We cannot directly deduce the tidal radius from our 
measurement. As the stellar density is still decreasing at 
the edge of the observed field (r ^ 110"), the King model 
can satisfactory fit for any tidal radius of r > 110". This 
difficulty is also expected from the fact that the EFF 
model, which does not take into account a tidal radius, 
fits the density profile equally well as the King model. 
We discuss other constrains for the tidal radius and the 
cluster size in the following section. 



8.2. Total mass of NGC 3603 
8.2.1. Size of NGC 3603 

In order to derive the total mass of the cluster by in- 
tegrating the radial mass density profile, the size of the 
cluster has to be known. Since fitting the King model 
to the density profile cannot realistically constrain the 
cluster's tidal radius, we adopt a cluster size from liter- 
ature. NPG02 give r = 150" ± 15" as the radius where 
the stellar density falls below 3 cr of the background vari- 
ation. They also give a tidal radius rt ^ 1300" from the 
best-fit King model, but as in our study, this radius is 
outside their ob served field, and accordingly very uncer- 
tain. Similarly, ISung fc Besselll (|2004D report a cluster 
radius of 2' based on the density profile of the bright 
stars identified in their optical and X-ray observations. 
They derive a tidal radius of rt ^ 900". Following these 
two studies, we adopt a cluster size of r ~ 150". 

In addition, we give upper and lower limits of the clus- 
ter size to properly estimate the total mass and half-mass 
radius of the cluster with the uncertainties. First, given 
the steady decrease of the density profile within our field, 
we can simply give a conservative lower limit of r = 110". 
The upper limit is derived from the cluster's tidal radius 
using the Galactic rotation curve and the measured den- 
sity profile in a self-consistent way. I n the study of the 
Pleiades cluster iPinfield et all (|1998f ) give an approxi- 
mate expression for the tidal radius of the cluster. 



GMc 



1/3 



- 1 46A^^/^ 
2{A-B)\ " ' 

where A4c is the mass of the cluster, and {A — B) is 
the difference of the Oort constants, which describes the 
differential rotation of the Galaxy in t he solar neighbor- 
hood, from lKerr fc Lvnden-BeTl (|l986f ). Since the Galac- 
tocentric distance of NGC 3603 is similar to that of the 
Sun, we adopt this relation for our estimate. Using equa- 
tion ([5]) and the total mass calculated from the mass 
density profile (which will be explained below), we can 
iteratively determine a self-consistent tidal radius. As a 
result we derive rt 1260" for the upper limit of the 
cluster size. 

8.2.2. Deprojection of the surface mass density distribution 

Integrating the best-fit King model, in which our esti- 
mate of the tidal radius rt — 1260" is set, up to the three 
radii r = 110", 150", and 1260", we derive the best esti- 
mate of the total mass and half-mass radius with error 
bars. Here we use a volume density distribution derived 
from the deprojection of the best-fit King model. The 
calculation is done in the following way: (1) deprojection 
of the King model into a volume density, (2) integration 
of the volume mass density up to the cluster radius, and 
(3) scaling the stellar mass from the mass range used for 
the mass density (0.5 — 2.5 Mq) to the full stellar mass 
range (assuming 0.1 — 100 Mq) using the derived IMF. 
In a spherical system, the surface mass density f{R) is 
related to the volume mass density m{r) by the Abel 
integral 

/(i?) = 2/J^S^. (6) 



(5) 



i?2 



Inverting this equation gives the volume mass density 
1 df{R) dR 



m{r) 



dR Vi?2 - r2 



(7) 
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where r and R are the spatial and the projected radii, 
respectively. The enclosed mass in a sphere of radius r 
is then obtained by integrating equation ([7]) as A^(r) = 
47r m{x)x^dx. 

First we derive the upper limits. Integration of the 
de-projected King model up to the tidal radius of 1260" 
gives a total mass of ~ 2520 Mq for the stars within 
the mass range of 0.5 — 2.5 Mq. Assuming a power law 
IMF with r = -0.74 over the mass range 0.1 - 100 Mq, 
we then obtain a total stellar mass of Aitot ^ 16,000 
Mq. The half-mass radius of NGC 3603 is rhm ~ 52" 
(~ 0.7 pc at 6 kpc). As we have measured the radial mass 
density profile for stars within 0.5 — 2.5 Mq, this half- 
mass radius applies only to stars in this mass range and 
does not take into account any mass segregation. For the 
best-guess radius of r = 150", we derive Mtot ^ 11800 
Mq and rhm ~ 30". For the lower limit of r = 110", we 
derived Mtot 10700 Mq and rhm ^ 25". 

To get a handle on the uncertainty of the derived to- 
tal mass, we analyze the potential effect due to system- 
atic errors in the assumed IMF power law index and the 
stellar mass range. As our goal is to estimate the total 
stellar mass outside the observed field of view r > 110", 
it may be more suitable to use the IMF of the outermost 
ISAAC fields (F = -0.85) instead of the average IMF. 
In this case we derive Mtot ~ 9700 Mq, i.e., an ~ 18% 
decrease. Next, if we assume a slightly wider mass range 
of 0.01 - 120 Mq instead of 0.1 - 100 Mq, we derive 
Mtot ^ 13500 Mq, an ~ 15% increase. Taking into ac- 
count the lower and upper limits for the cluster radius 
and the errors from the IMF slope and the stellar mass 
range, we estimate the total mass of NGC 3603 to be 
1.0 - 1.6 X lO'^ Mq, and the half-mass radius to be about 
25" - 50" (0.7 - 1.5 pc at d = 6 kpc). 

Using the core radius of ~ 5" we also derive the central 
mass density of ~ 6 x 10'' Mq pc~^. It is worth men- 
tioning that the core density measured directly from the 
surface density profile (without deprojection) would give 
a factor of ~ 2 larger value. 

8.3. Dynamical evolutionary status 

Here we estimate the relaxation time of the cluster. As 
the relaxation time depends on the stellar density, it is a 
function of radial distance from the cluster center, and it 
varies by several orders of magnitude across the cluster. 
While the inner parts relax most quickly, the outer parts 
take much longer to relax. The relaxation time also de- 
pends on the stellar mass. Higher-mass stars relax more 
quickly than lower-mass stars. In our study we restrict 
ourselves to the average relaxation time of the cluster. 
Following i Binnev fc Tremaind (|1987f ). the median relax- 
ation time trei of a stellar population is given by 

_ 6.5 X 108 f M \ / IMq \ / reh \ 
'''' ~ ln(0.47V) \10^ Mq J J\lpcJ 

where M is the total mass within the characteristic ra- 
dius rch, rn^ is the characteristic stellar mass, and N 
is a total number of stars in the cluster. Usually the 
half-mass r adius rhm is taken as the characteristic ra- 
dius (e.g., Portegies Z wart et al.ll2002l ). The according 
half-mass relaxation time is then a good estimate for the 
dynamical timescale of the whole cluster. 



Here we estimate the relaxation time for the two ex- 
treme cases with upper and lower limits of the total mass 
and half-mass radius. For the upper limit half-mass ra- 
dius rhm = 50", a total mass of ~ 8000 Mq within the 
radius, a characteristic stellar mass of 1 Mq, and a total 
of ~ 16,000 stars, we get a half-mass relaxation time of 
~ 39 Myr. As a lower limit, we derive ~ 11 Myr. The un- 
certainty in the number of stars could change the relax- 
ation time by ±5 — 10%. In any case it is safe to say that 
the dynamical timescale for stars with masses around 1 
Mq is of order 10 million years. This is about an order of 
magnitude larger than the age of the intermediate- and 
low- mass stars (< 2 Mq) in the cluster, which is less than 
~ 1 Myr. Thus, solar-mass stars are still young enough 
such that the dynamical evolution has not yet changed 
substantially their initial kinetic energy. 

As the dynamical time scale is inversely proportional 
to the stellar mass, the high-mass population is expected 
to have a much shorter relaxation time. The smaller half- 
mass radius of the high-mass stars further decreases the 
relaxation time proportionally to the power of 1.5 (see 
eq. [8]). For example, if we simply apply the lower limit 
of rhm = 25", which is derived based on the 0.5 — 2.5 Mq 
stars, we obtain a rough estimate of the relaxation time 
of about 1 Myr for stars typically of 10 Mq. We there- 
fore conclude that very high-mass stars in the central 
cluster in NGC 3603 have a dynamical timescale com- 
parable to the cluster age. Thus, we cannot conclude if 
the observed mass segregation in the high-mass popula- 
tion is dynamical or primordial. In fact it might be due 
to the mixture of both effects, since the position of mas- 
sive stars in young clusters generally reflec ts the cluster's 
initial conditions (|Bonnell fc Davieall998l ). A similar dy- 
namical state is seen in other young stellar populations. 
For example, the young ONC (< 1 Myr) has been re- 
ported to show a mass segregation, in which higher mass 
stars are preferentially located in the cluster center due 
to primordial effects. But as in NGC 3603, there is no 
evidence of mass segreg ation for stars below 1 — 2 Mq 
(|Hillenbrand et al.lll998l ). 

9. SYSTEMATIC UNCERTAINTIES OF THE IMF 
DETERMINATION 

There are some systematic uncertainties that could po- 
tentially affect the characteristics of the derived IMF. 
A first class of such uncertainties could arise from the 
age, distance, and foreground extinction of the cluster. 
The selection of stellar evolutionary model among vari- 
ous available choices is also subject to an uncertainty, and 
once a set of evolutionary model is chosen, the metallic- 
ity selection is also a potential uncertainty. From a tech- 
nical point of view, our incompleteness correction and 
the selection of cluster stars are potential uncertainties. 
Moreover, the stellar mass outside the observed field of 
view is a source of uncertainty. As for the intrinsic prop- 
erties of the cluster, the presence of variable extinction 
and unresolved binary/multiple systems could also affect 
the derived IMF characteristics. We thus investigate how 
these uncertainties could vary the index of the power law 
slope by scrutinizing them one by one. 

For the sake of simplicity, we use several techniques de- 
pending on the parameters. To examine the uncertainties 
due to age, distance, foreground extinction, metallicity, 
evolutionary model, and the stellar mass outside the ob- 
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served field, we construct a simplified two-mass bin IMF 
in each case. Here we count only the stars in two mass 
ranges, 1.7 — 2.6 Mq and 8 — 19 Mq as low- and high- 
mass bins, respectively. Another simplification is that we 
use only the J-band magnitude for the luminosity-mass 
conversion, in contrast of using all three JHKs bands 
for the best estimate of the IMF in § [71 For investigating 
the rest of the uncertainties, we follow the technique used 
in the determination of the best IMF. As some of the er- 
rors have asymmetric distributions, we give the median 
and the upper and lower limits in each measurement. In 
the end we combine all resulting errors to obtain the un- 
certainty in the best estimate of the power law of the 
IMF. 

9.1. Age 

As presented in §[521 we have used the 2.5 Myr LSOl 
isochrone for the MS population and the 0.7 Myr PS99 
isochrone for the higher mass PMS population (> 1.2 
Mq). To safely cover the uncertainties in the current age 
estimate of < 3 Myr, we study the effect of a selection 
among the 0.2, 1, 2.5 and 5 Myr LSOl isochrones for 
the high-mass bin, and the 0.3, 0.7, 1, and 1.5 Myr PS99 
isochrones for the low-mass bin. By counting the number 
of stars in the two mass bins, we compute the IMF power 
law indices for all possible combinations of the MS and 
PMS isochrones. We summarize the results in Tab. [3l 
There are several points to highlight. First, there is a 
correlation between the age and the resulting power law 
index. If the age of the MS isochrone is fixed to 2.5 Myr, 
a younger age of the PMS isochrone results in a shallower 
IMF slope. In contrast, for a fixed PMS isochrone age, 
a younger age of the MS isochrone results in a steeper 
slope of the IMF. This outcome can be explained by re- 
writing the definition for the IMF using a M.-L 
relation 

^{\ogM)^^{M)M(^^ , (9) 

where the last term is the derivative of the M-L rela- 
tion. In the M-L diagram (log — Mj) of the various 
isochrones, a younger age shows a steeper profile. Hence, 
a given mass range covers a narrower magnitude range for 
younger isochrones, resulting in smaller number counts. 
This leads to the observed correlation. Another point is 
that the biggest change of the index originates from the 
selection of the PMS isochrone. This is because the M.-L 
relation of PMS stars is more sensitive to the age varia- 
tions than that of MS stars. We derive the most extreme 
change for the youngest 0.3 Myr PMS isochrone. 

We note that the IMF power law index derived from 
the two mass bins is only little different from fitting the 
complete IMF. For our best age estimate with 2.5 Myr 
MS and 0.7 Myr PMS isochrones, the index from the sim- 
ple two-mass bin analysis gives F — —0.82 in comparison 
to r = —0.74 from the complete fit. 

As a result, we find that the age uncertainty leads to 
a systematic uncertainty in the IMF power law index of 
AF = -1-0.60 and —0.41. We note that the systematic 
error in the above analysis provides the upper and lower 
limits, but a realistic uncertainty is expected to be much 
smaller. We are aware that the 1.5 Myr PMS isochrone 
clearly fails to fit the MS turnon point and that the 0.3 
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Variation of the power law index due to age 
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of isochrone 
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tligh-mass bin 
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0.2 




0.3 
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1.0 




0.3 


-0.34 
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-0.30 


5.0 
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-0.28 


0.2 
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-0.87 


1.0 




0.7 


-0.86 
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-0.82 


5.0 




0.7 


-0.80 


0.2 




1.0 


-1.10 


1.0 




1.0 


-1.09 


2.5 




1.0 


-1.05 


5.0 




1.0 


-1.03 


0.2 




1.5 


-1.28 


1.0 




1.5 


-1.27 


2.5 




1.5 


-1.23 


5.0 




1.5 


-1.21 



Myr PMS isochrone fails to reproduce the bulk of the 
population in the PMS-MS transition region. 

9.2. Distance 

Earlier studies have measured the distance to be within 
6 — 8 kpc, and d ~ 6 kpc yields a reasonable fit in our 
isochrone fitting in the CMD. We therefore study dis- 
tances of 5, 6, 7, and 8 kpc to see the resulting change of 
the IMF slope. Note that we do not adjust the color cut 
for the rejection of field stars but only change the dis- 
tance in the apparent magnitude-mass conversion. The 
resulting power law indices are F = —0.85, —0.84, —0.86, 
and —0.85 for d = 5 — 8 kpc, respectively, with a neg- 
ligible variation of AF = ±0.01. The small variation 
of the index is expected from the fact that a change in 
distance simultaneously shifts both low- and high-mass 
magnitude ranges, and thus does not yield any significant 
differences in the ratio of the high- to low-mass counts. 

9.3. Foreground extinction 

As described in 15.31 our best estimate of the interstel- 
lar extinction towards NGC 3603 is Av — 4.5 ±0.5 mag. 
Considering the observed radial increase of Av4y ~ 2.0 
mag from the cluster center to the outermost ISAAC 
fields, we apply three different values Ay = 4.0,4.5, 
and 6.5 mag to derive the uncertainty in the IMF slope. 
Again we do not adjust the color cut for the rejection of 
field stars. The resulting power law indices of the IMF 

are F 0.86,-0.83,-0.86 for Ay = 4.0,4.5,6.5 mag, 

respectively. The slope does not show any trend and is 
fairly insensitive to the adopted foreground extinction 
with an error of only AF = ^q'qi- 

9.4. Theoretical model dependence 

Here we adopt several currently available MS and PMS 
evolutionary models to derive the number counts for the 
low- and high-mass bins and calculate the change in the 
power law index of the resulting IMF. 

For the 8—19 Mq high-mass bin, in addition to the 
2.5 Myr MS isochrone in the Geneva model (LSOl), we 
apply an isochrone of the Padova evolut ionary model 
computed with the 2MASS filter systern (iBertelli et "aTl 
11994 iGirardi et al.ll2000t iBonatto eral]|2004f ). We use 
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TABLE 4 

Variation of the power law index due to the selection of 



STELLAR evolutionary MODELS 


Isochronc model 


r 


High-mass bin (MS) Low-mass bin (PMS) 


LSOl (Geneva) 2.5 Myr PS99 0.7 Myr 
SDFOO 1.5 Myr 
y2 1.5 Myr 


-0.82 
-0.75 
-0.95 


Padova 4 Myr PS99 0.7 Myr 
SDFOO 1.5 Myr 
y2 1.5 Myr 


-0.83 
-0.76 
-0.95 



the youngest available 4 Myr isochrone here. For the 
1.7 — 2.6 Mq low-mass bin, in additio n to the PS99 set , 
we apply a PMS model pubhshed bv ISiess et al.l ()200ClL 
hereafter SDFOO), and the isochrones based on the re - 
vised Yale models ()Yi et al.ll200lUDemaraue et al.|[200l . 
As discussed in § 15. 2( it is well known that the various 
available PMS evolutionary models result in considerable 
discrepancies in the age estimate of a star cluster, in par- 
ticular, for a very young age (less than a few megayears). 
Although we have applied the 0.7 Myr PS99 isochrone for 
the PMS stars, for this analysis, it is necessary to individ- 
ually select the best-fit age for each PMS model instead of 
simply applying 0.7 Myr. We select the 1.5 Myr SDFOO 
isochrone and the 1.5 Myr isochrone for the PMS. 
The resulting IMF power law indices are summarized in 
Tab. m showing a slightly asymmetric distribution with 
AF 



4-0.07 
-0.13- 



9.5. Metallicity 

We have so far adopted a solar metallicity Z — 0.02 
in the isochrone selections. This is based on the fact 
that the Galactocentric distance of NGC 3603 is sim- 
ilar to that of the Sun, so that the radial gradient of 
the metallicity in the Milky Way does not n eed to be 
taken into account (e.g. iRudolph et al. I [20061 ). To esti- 
mate the potential effect of the metallicity selection, we 
test a half-solar, solar, and twice-solar metallicity. Since 
no isochrones with different metallicities are available in 
the PS99 set, we employ the 1.5 Myr SDFOO PMS models 
with z — 0.01, 0.02, and 0.04 for the low-mass stars. For 
the high-mass stars we use 2.5 Myr LSOl MS isochrones 
with z = 0.008, 0.02, and 0.04. The resulting power 
law indices are F = —0.88 (half-solar), —0.75 (solar) and 
—0.76 (twice-solar). The application of the solar and 
twice-solar metallicity yield almost the same value, and 
a major difference is seen for the half-solar metallicity, 
which steepens the IMF by AF ^ -0.13. Note that, al- 
though we confirmed that there were clear correlations 
between the metallicity and the number counts in the 
both low- and high-mass bins, no trend eventually was 
seen in the power law index. Although we do not find any 
trend in our analysis, a possible dependence of the IMF 
of stellar systems on the met allicity has been discussed 
in several recent studies (e.g.. iSoUima e t al.. 2007.) . 

9.6. Individual extinction 

We have so far applied a uniform foreground extinction 
of Ay — 4.5 mag in the conversion from the observed 
luminosities to stellar masses. However, potential local 
variations from circumstellar gas and dust, and larger 



scale variations on the size scale of the cluster could af- 
fect the derived stellar masses. In fact we find a radial 
variation of AAy ~ 2.0 mag from the cluster center to- 
wards the outer regions. To derive the impact of this 
uncertainty, we construct the IMF of NGC 3603 with 
and without the treatment of a variable extinction. For 
this we shift the stars along the reddening direction in 
the JHKs magnitudes space as close as possible to the 
isochrone, rather than picking the closest points on the 
isochrone for a fixed average extinction. 

As a result, we find that the treatment of the vari- 
able extinction flattens the IMF slope by AF ^ -1-0.15. 
We also find that the change of the IMF slope is mostly 
caused by the almost parallel alignment of the reddening 
direction and the isochrone of the PMS-MS transition re- 
gion. Because of this alignment, sources in the scattered 
PMS-MS region are dereddened to the intermediate-mass 
MS, which causes an artificial bump at around 5 Mq and 
results in a slightly flatter IMF slope within 0.4 — 20 Mq. 
As the change of the power law index mostly arises from 
this local feature, the real impact of the presence of the 
individual extinction on the IMF slope is expected to be 
even smaller. 

9.7. Incompleteness correction 

Here we evaluate the potential influence of our 
empirical-based incompleteness correction on the IMF by 
considering two extreme cases: one is without any correc- 
tion, and the other is with an over-correction. The first 
case has already been analyzed in Figure 1121 the power 
law index being F = —0.68. This means that without the 
incompleteness correction the IMF would come out shal- 
lower by AF ^ -1-0.06. For the second case, we prepared 
a larger correction option in the computation of the cor- 
rection factor (see detailed explanation in § [4|) . Using the 
same technique as for the best IMF (F = —0.74), but ap- 
plying a larger correction factor, we find the IMF showing 
a power law index of F ~ —0.87. Thus, the larger cor- 
rection steepens the IMF by AF ^ —0.13. Here we note 
that a correction by simply a factor of 2 of the normal 
correction in each mass bin steepens the IMF power law 
index by AF 0.04. 

In summary we derive a potential error of AF =1^^ 
in our incompleteness correction. As mentioned in § 17.31 
the slight differenc e of the resulting IMF between our 
study and that by IStolte et al, (2006) would be partly 
explained by the difference in the correction rates be- 
tween two studies. 

9.8. Cluster membership 

In order to deal with line-of-sight field stars in the ob- 
served fields, we have so far applied the color cut method 
instead of a statistical approach based on control field 
measurements. In the color cut we exclude blue stars 
in the (J — Ks,J) CMD because those stars can most 
likely be regarded as non-cluster members. As mentioned 
in § 16.11 a potential source of uncertainty in this applica- 
tion is the fact that the region redward of the color cut 
could still contain some field stars whose J — Ks colors 
are either similar to those of cluster stars, or are redder 
than cluster stars, e.g., red giants^. 

* We note that extragalactic contaminations are expected to 
be negUgible because of the relatively small field of view. Quan- 
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In fact, we verify that the "blueward only" color cut 
is a reasonable treatment based on the analysis of LF 
in § O adopting the field star KLF in NPG02 which is 
derived from a control field and is in agreement with what 
is expected from Galactic models. Although we arc not 
able to make any statistical estimate of field stars from 
our own data because of the lack of a control field in our 
observations, we are still able to test the potential effect 
of red stars by simply cutting them out ("red cut"). 

We determine the red cut criteria as follows. For 
MS stars down to mj < 14 mag, stars located within 
A( J — Ks) = +0.35 mag redward from the combined 
isochrone are regarded as cluster stars, and stars redder 
than the boundary are excluded from the IMF determi- 
nation. For PMS stars with mj > 15 mag, we apply 
A( J — Kg) = +0.50 mag considering larger photometry 
errors for fainter stars. As for the region of mj — 14—15 
mag, the line is drawn by directly connecting the curves 
in the upper and lower magnitude ranges so as to deal 
with the turnover feature of the isochrone in the PMS- 
MS transition region. In addition, since we measure the 
radial increase of the interstellar extinction across the 
cluster (as shown in Figure El), we separate inner and 
outer regions with r < 30" and r = 30" — 110" and ap- 
ply an additional shift of +0.25 mag for the outer region. 
Thus, we use A( J — Kg) = +0.60 and +0.75 mag for 
MS and PMS stars, respectively. Here we note that, as 
the blueward color cut, the shift is applied only in the 
J — Ks direction. 

The red cut is illustrated in the (J — Kg, J) CMDs in 
Figure [T51 As a result, 810 stars out of 7514 stars, which 
are used for the best estimate of the IMF, are classified as 
red stars and excluded. Thus, 6704 stars are used for the 
IMF calculation, and the resulting IMF shows the power 
law index of F = —0.76, being fairly similar to that of 
the best IMF. Thus we conclude that the potential error 
of the IMF power law index due to the presence of red 
stars is AF 0.02. 

We note that our criteria of the red cut i s simi lar to 
those in the previous study bv IStolte et al.l (l2006l ). and 
the only difference is that we apply two different thresh- 
olds depending on the radial distance. 

9.9. Unresolved binary systems 

Unresolved binaries (or higher-order multiple sys- 
tems) have been identified as a n important uncer- 
tainty in the IMF determ ination (jKroupal 1200 iL l2002t 
iMalkov fc Zinneckeill200lh . Any unresolved source com- 
posed of two or more stars mimics a slightly brighter sin- 
gle star in the photometry, resulting in an underestimate 
(overestimate) of the low-mass (high-mass) population. 
A mass function derived from a population with unre- 
solved systems, the so-called system mass function, thus 
shows a slightly larger power law index than the actual 
stellar mass function, and an appropriate treatment of 
the unresolved sources will steepen the derived IMF. 

In this section we first present our estimate of the bi- 
nary fraction based on the analysis of the (J — A'^, J) 
CMD. In a second step we then model the impact of 

titatively, integrated galaxy number densities are expected to be 
approx. 10* galaxies deg ~'^ up to ~ 20 mag in t he H and K 
bands (e.g. Dior govski et al~lll995l : lYan et al.lll998l ). correspond- 
ing to only a few galaxy counts in our field of r < 110" 



the unresolved binaries on the IMF by applying the well 
known binary frequency from the ONC. 

9.9.1. Unresolved binary fraction from the observed CMD 

Several studies have already reported the potential bias 
from unresolved binaries on the IMF of Galact ic young 
star clusters. For example, iFiger et all (|1999( l derived 
a significantly fiat IMF with F ~ —0.65 for the Arches 
cluster, for which a binary fraction of unity would steepen 
the IMF by AF 0.3. 

Stolte et all (|2004 l2006f) report the presence of a sec- 
ondary sequence in NGC 3603, which shows up in their 
{J -Ks, J) CMD towards the bright red side of the MS 
above the PMS-MS transition region. They interpret 
this sequence as equal-mass binary stars, formed together 
with the other PMS stars in a single star formation epoch 
~ 1 Myr ago. They find no IR excess emission for the 
stars of this secondary sequence in their L-band analysis 
but rather find an anti-correlation between those sources 
and H-a emission stars. This suggests that the secondary 
sequence is not caused by the contribution from circum- 
stellar disks, but by unresolved binaries. Counting stars 
on the secondary sequence, they derive a fraction of un- 
resolved binaries of 30% in the central cluster. From the 
typical offset of the secondary sequence of A J ~ —0.75 
mag, they argue that most of the candidates are likely 
to be close to equal-mass binary systems. Using the 
unresolved binary fraction of 30% and assuming equal- 
mass components, the authors apply a binary correction 
in the IMF determination, steepening the IMF slope by 
AF 0.06. 

In contrast, we do not find a distinct secondary se- 
quence, but merely a scattered distribution towards the 
red side of the MS and to the top of the PMS-MS transi- 
tion region. Therefore we think it is not feasible to confi- 
dently distinguish potential unresolved binaries (which 
can have non-uniform mass-ratio) from reddened MS 
stars or younger PMS stars. The scatter can be partly 
due to an intrinsic age spread in NGC 3603 as discussed 
in our age estimate in § 15.21 

While it is difficult to reliably estimate the fraction of 
unresolved binaries from our data, we can give a rough 
estimate by counting sources located in the upper part of 
the PMS-MS transition in the {J-Ks, J) CMD. The se- 
lection criterion for potential binaries is illustrated in the 
CMDs in Figure [1^1 The location of the binary band is 
defined based on the following idea: in case of an equal- 
mass system, the real magnitude of the two stars is 0.75 
mag larger (i.e., fainter) than the system magnitude. As- 
suming a ±0.3 mag error in the J magnitude to define 
the width of the band, the limits of the binary band 
are set by shifting the isochrone by —1.05 and —0.45 
mag. Among the total of 204 sources at r < 13" within 
nij ~ 13—16 mag, we find 27 stars within the binary 
band, corresponding to a binary fraction of ~ 13% (Fig- 
ure [11] left). The same experiment for the outer ISAAC 
field of 13" < r < 110" yields a similar value of ^ 17% as 
96 binary candidates are found among 577 sources in the 
magnitude range. However, since this measurement is 
based on the assumption of equal-mass components and 
the typical uncertainty of ±0.3 mag, the true number is 
subject to large uncertainties. For example, if we adopt 
a magnitude range of Amj = 0.75 ± 0.4 mag to define 
the binary band, the binary fraction increases to ^ 21% 
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Fig. 15. — Selection of cluster stars by the red cut in addition to the blue cut for the inner field of r < 30" (left) and outer field of 
r = 30" — 110" (right). The gray solid curve redwards of the isochrone is the boundary of t he re d cut. Red stars regarded as non-cluster 
members and thus removed for the IMF determination are shown in gray triangles. As Figure [TO] blue stars based on the original color cut 
are shown in gray squares. 



for the NACO field and - 22% for the ISAAC field. 

Here we note that the potential unresolved binaries in 
our data also include multiple sources from projection 
effects, i.e., optical binaries. Using the spatial resolution 
of - 490 AU in the NACO field and - 2200 AU in the 
ISAAC field, we derive a likelihood for optical binaries 
of about 6 and 9% for the NACO and ISAAC data, re- 
spectively. 

9.9.2. Systematic error of the IMF from unresolved 
binaries 

As outlined above, we cannot reliably identify unre- 
solved binaries in our data in the CMD analysis. To 
quantify the systematic error of the IMF from such un- 
resolved binaries, we thus choose to simulate the effect 
adopting the measured binary distribution in the ONC 
considering the similar age and high-mass stellar content 
of the ONC and NGC 3603. 

Binary fraction of high- and low-mass stars in the ONC. — 
The binary fraction in the ONC including the Trapezium 
cluster is among the best studied. Note that, although 
some binary/multiple surveys mentioned below differen- 
tiate between binary and higher-order systems, we do 
not follow thi s distinction but sim ply use the term bi- 
nary fraction. iProsser et al.l (|1994D . for example, derived 
a binary fraction of ~ 12% for a projected linear sep- 
aration of 26 — 4 40 AU in the Trapezium cluster, and 
iPetr et al.l (fT998h derived ~ 9.5 ± 4.5% for 63 - 225 AU 
for stars with masses down to ^ 0.04 Mq. Combining 
their n ew observations with the results from iPetr et al.l 
(|1998[ ). iKohler et all ()2006D derive a binary fraction of 
~ 14 ± 7% and ~ 23 ± 10% for stars with M > 2 Mq 
in the periphery and the central core, respectively, and 
- 4.8±1.8% and - 3.6±3.2% for stars with X = 0.1-2 
Mq for a projected separation of 60—500 AU. The higher 



binary fraction for high-mass stars in the ONC is co nsis- 
tent with earlier results from Prcibis ch et al.l (|1999t). 

Her e we adopt the binary fractions from iKohler et ahl 
(|2006f ). To simplify our analysis, we use their typical 
values of ~ 5% for low-mass stars (0.1—2 Mq) and ~ 15% 
for intermediate- to high-mass stars (> 2 Mq). To also 
cover the various other estimates of the binary fraction in 
young stellar clusters and to estimate the upper limit of 
the potential impact of unresolved binaries on the IMF 
determination, we also analyze the case for a 2 times 
larger binary fraction, i.e., ~ 10% for 0.1 — 2 Mq and 
- 30% for > 2 Mq. 

Binary fraction as a function of separation. — The an- 
gular resolution of our NACO and ISAAC data cor- 
respond to minimum detectable binary separations of 
^ 490 and ^ 2200 AU, respecti vely. In their binary 
analysis of nearby solar-type stars. iDuauennov fc Mavon 
(1991) show that the orbital period f ollows a lognormal 
distribution. Using this distribution iKohled ()2001l ) de- 
rives the distribution for the projected separation, which 
follows the form f{x) oc exp[— i(a; — x)'^/a'^], where x 
is the logarithm of the projected separation in AU with 
X = 1.44 and a — 1.55 (see Figure 6 in the reference). 
Integrating this distribution up to the minimum sepa- 
rations, we can estimate that 79% of all binaries are 
unresolved in the NACO data, and ^ 89% are unresolved 
in the ISAAC data. Considering the rather larger uncer- 
tainty in the adopted binary fractions, we ignore this 
difference in the following experiment. 

Binary fraction as a function of mass ratio. — The mass 
ratio {q — m2/mi) distribution has shown various forms 
in var ious types of stellar populations. iReid fc GizisI 
(|1997f ) find evidence for preferential formation of nearly 
equal- mass binary systems, i.e., a peak at q ^ 1 among 
nearby M dwarfs. In contrast, IDuauennov fc Mavorl 
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Fig. 16.— Unresolved binary candidates in the (J - Ks, J) CMDs for stars from the NACO field of r < 13" {left) and the ISAAC field 
of 13" < r < 110" {right). The analysis is restricted to a range of mj = 13 — 16 mag. The binary band is indicated in gray solid curves. In 
total 186 cluster stars {large dots) plus 18 non-cluster stars {plus signs) are identified in the magnitude range, and 27 sources are located 
within the binary band {points and squares), i.e., a binary fraction of ~ 13%. In the outer field, 96 binary candidates among 577 stars are 
detected, resulting in ~ 17%. 



(|1991| ) derive a distributi on that peaks near q = .2 — 0.3 
for solar-type field stars. iFischer fc Marcvl ()1992l ) derive 
a peak at q = 0.4 — 0.6 and 0. 6 — 0.8, respectively, fo r the 
field G and M dwarf binaries. iDuchene et al.l (119991 ) find 
a distribution peaking at q = 0.5—0.75 for binaries with a 
primary mass of < 1 M© in the young star cluster IC 348. 
For some clusters the mass ratio distribution can be fitted 
by an inverse power law f(q ) oc q~^ , for example, with 
r = 0.33 - 0.5 in Sco 0B 2 (|Shatskv fc TokovininI [200l 
iKouwenhoven et all 120051) a nd F = 1.3 in the Hyades 
cluster ( Patience et al. I TT998I) . 

Thus, it is not straightforward to assume a unique mass 
ratio distribution; we hence examine two cases. In the 
first model, which favors equal-mass binaries, 75% of the 
unresolved binary systems have a mass ratio of q — 1 
and 25% have q = 0.4. In the second model, which favors 
non-equal-mass binaries, 25% of the population have q = 
1 and 75% have q = 0.4. 

Simulation of the unresolved binary correction. — Using 
the adopted binary fractions and the two cases of the 
mass ratio distribution mentioned above, we study the 
potential influence of the unresolved binaries on our IMF 
determination by means of a simple simulation. We per- 
form the simulation as follows. First we randomly gen- 
erate a set of binaries for each mass bin according to the 
observed number counts and the assumed binary frac- 
tions. We then split those binary systems according to 
the mass ratio distribution, subtract the system from the 
original mass bin, and then add the split masses into the 
corresponding mass bins. Finally we compute the power 
law index of this simulated IMF within the mass range of 
0.4 — 20 Mq. As a result, the power law index of the IMF 
becomes slightly steeper by AF ~ —0.02. Even assuming 
twice the binary fraction of the ONC, the IMF steepens 
by only AF ~ —0.04. The two cases of the mass ratios 
yield almost identical power law indices. We conclude 
that the systematic error in the IMF power law index 
from the presence of unresolved binaries is AF < —0.04. 

9.10. The stellar mass outside the observed field of 



In the analysis of the dynamical state of the cluster in 
§ m we have outlined that the high-mass stars in NGC 
3603 appear to be strongly concentrated towards the 
cluster center. In contract, the intermediate- and low- 
mass stars in the cluster appear neither dynamically nor 
primordially mass-segregated at this stage and extend 
beyond the observed field of r < 110". Therefore, we 
need to quantify how much these intermediate- and low- 
mass stars outside the field of view can change the IMF. 

For that we use the de-projected mass density distribu- 
tion from § 18.21 Integrating the volume mass density of 
the stars within the mass range of 0.5 — 2.5 Mq up to the 
observation limit r = 110", we find that our observations 
cover more than 90% of these stars if the cluster radius 
is r = 150", and still ~ 67% when using the upper limit 
of the cluster size, i.e., the tidal radius of vt = 1260". 
Qualitatively we conclude that our data indeed covers the 
majority of the cluster mass, and that the observed IMF 
is thus representative for the whole cluster. To quantify 
how much the intermediate- and low-mass stars outside 
our field of view can change the IMF power law index, 
we simply add the maximum missing mass of 33% to the 
number counts of the low-mass bin (1.7 — 2.6 Mq) in 
the simplified two mass bins IMF. As a result, the power 
law index decreases from F = —0.82 to —0.98, that is, a 
steepening of the IMF with AF -0.16. 

9.11. Error combination 

Tab. [5] summarizes the individual systematic errors 
from l9.1l - [97T0l We combine them into a single system- 
atic error to obtain the final power law index in our IMF 
determination. To c ombine those a symmetric errors, we 
apply the method bv lBarlowl (|2003D ^ The program con- 
structs, for each set of positive and negative errors, an 
asymmetric distribution consisting of two half-Gaussian 
distributions and then convolves these distributions to 
derive a combined asymmetric error. We note that in this 
combination we also add the fitting error ±0.02 in the 

^ This method is implemented in the program from 
|http: //www. slac. Stanford. edu~barlow/java/statisticsl.html| 
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TABLE 5 

Summary of the error estimates and the combination 



PfircimGtGrs 






Age 


fin 




Distcirice 


±0.01 




Foreground extinction 


+0.02, 


-0.01 


Evolutionary model 


+0.07, 


-0.13 


Metallicity 


-0.13 




Individual extinction 


+0.15 




Incompleteness correction 


+0.06, 


-0.13 


Cluster membership 


-0.02 




Unresolved binary 


-0.04 




Stellar mass outside the FOV 


-0.16 




Combined 


+0.62, 


-0.47 



power law index calculation. The resulting combined er- 
ror is Ar = !to'47- For comparison, the simple quadratic 
sum of the positive and negative errors yields fairly sim- 
ilar values of -1-0.63 and —0.50, respectively. 

In summary, we conclude that the IMF power law in- 
dex of NGC 3603 for the mass range 0.4 - 20 Mq in- 
cluding all systematic errors is F = — 0.74^Q'4y. Consid- 
ering a Gaussian probability distribution with the above 
1 a errors, we calculate the probability that the IMF of 
NGC 3603 is as steep as the Salpeter IMF, resulting in 
Pr(r < -1.35) ~ 0.1, i.e., about 10%. This leads us to 
interpret that the IMF of NGC 3603 is distinguishably 
flatter than the Salpeter- like IMF. 

10. DISCUSSION 

We find that the IMF of NGC 3603 is weU described by 
a single power law with F = -O.TAto 'ij within 0.4 - 20 
Mq. We also find that the IMF steepens towards outer 
region with the power law index from F ^ —0.3 to —0.8. 
The steepening mainly occurs in the inner r < 30". The 
average IMF is slightly flatter but still consistent with 
previous studies (see discussions in § I7.3P . Here we dis- 
cuss the resulting IMF of NGC 3603 with a particular 
focus on the question: is the IMF of NGC 3603 an ex- 
ample for top-heavy IMF in starbursts? 

10.1. IMF in various stellar populations 

There is growing evidence that the IMF varies between 
different stellar populations. Here we summarize some 
recent IMF studies of young stellar populations. 

Several stellar clusters have been suggested to have a 
somewhat flat IMF. For example, the Arches cluster near 
the Galactic Center (GC) shows a comparatively shal- 
low IMF with r = —0.6 to -1. 1 for the high-mass star s 
(IFieer et all [19991 : IStolte et al.l 120051: iKim et al.ll200l . 
The Quintuplet cluster - another massive young clus- 
ter near the GC - has also been suggested to have a 
mass distribution potentially flatter than the field IMF 
(IFieer et al.lll999D . 

Other stellar clusters show standard Salpeter-like 
IMFs, but there are substantial differences between var- 
ious studies. For example, the IMF of the massive star 
cluster R136 in 30 Dor has been report ed to have a 
powe r law index of F ^ ~\ -^ *o —1.6 (jBrandl et al.l 
119961 iMassev k Hunteii I1998D. and even subs tantiallv 
flatter for solar-mass stars (jSirianni et al.ll200Cil ). There 
are many other examples such as the double cluster 



h and y Persi with power law indices of F ~ —1.3 
(|Slesnick et al.l[200l . and NG C 1960 (M36) and NGC 
2194 with F 1.2 to -1.3 (jSanner et al.ll2000D . Re- 
cently, the IMF of the Galactic massive c luster Wester- 
lund 2 was measured to fit with F ~ —1.2 (lAscenso et al.l 
|2007[) . Yet other clusters have IMFs slightly steeper, 

for example, NGC 242 2 with F 2.0 (0.9 - 2.5 Mq) 

(jPrisinzano et al."2003'), NGC 3576 with F = -1.62 (> 3 
Mq) (Figuercdo et al. 2002) , and the Sco 0B2 associa- 
tion with r = -1.9 ()Brownlll998[ ). 

From studies of nearby star- forming regions, it has 
been found that the variations in the IMF power law in- 
dex and in the characteristic masses become even larger 
towards very low stellar masses. Examples of objects 
which show su ch variations are the Trapezi um cluster 
(iMuench et al.l l2002') , the Pleiades cluster (Bouvi er et al. 
1998), the Taurus star-forming region (Briccno et al. 
2002), the embedded cluster IC 348 (Preibisch et al. 
2003; Luhman et al. 200|), and the NGC 1333 molecular 
cloud (iWilking et al.ll200 4). 

Thus, we have so far known that there are moderate 
variations in the IMF characteristics in intermediate- to 
high-mass stars among young stellar clusters and sub- 
stantial variations at very low- and substellar masses. 

10.2. Universal or variable IMF? 

The question arises, is the observed IMF variation a 
true deviation from the Salpeter IMF? Or can the varia- 
tion be explained by systematic uncertainties in the mea- 
surements? 

[Kroupa i2001f ). for example, argues that the so- far ob- 
served IMF variations are caused by combined effects of 
observati onal, th e oretic al , and statistical uncertainties. 
Prior to iKroupal ()2001f ). iScalol (|1998l ) summarized the 
state of the IMF research by compiling the IMFs of many 
clusters and associations. Although no systematic trend 
is seen, there is a substantial scatter of at least unity 
(±0.5) in the F versus stellar mass plot, indicative of IMF 
variations. The large uncertainties prevent the author 
from giving any conclusive statement about a varying 
IMF. As an average IMF, the author proposes a three- 
segment power law IMF with the index r = — 0.2±0.3 
for the mass range of 0.1 — 1 Mq, F = —1.7 ± 0.5 for 
1 - 10 Mq, and F = -1.3 ± 0.5 for 10 - 100 Mq. Here 
t he IMF i s steep est in the intermediate-mass range. 

IKroupal (|200l| ) defines the Galactic field IMF adding 
the results from local star cou nts to the c ompilation 
of the MF power law indices in IScalol ()1998l ) . Assum- 
ing the defined average IMF, the author investigates 
the uncertainty inherent in any observational IMF es- 
timate through A^-body model calculations. The au- 
thor concludes that no true variations can be detected 
within the fundamental limit of uncertainties of a uni- 
versal IMF described by a segmented power law with 
index of F = 0.7 ± 0.7 for 0.01 - 0.08 M©, F = -0.8 ± 0.5 
for 0.08 - 0.5 Mq, T ^ -1.7 ± 0.3 for 0.5 - 1 Mq, and 
F = -1.3 ± 0.7 fo r > 1 M q. 

Recently IScalol ()2005l ) again summarized the current 
IMF estimates from field star counts and from observa- 
tions of open clusters. The author found that the vari- 
ation of the IMF power law index among the field star 
studies within 1 — 15 Mq is not negligible and that the 
variation among the cluster IMFs is considerable. These 
variations might be accommodated by the combined ef- 
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fects of various study-specific uncertainties among stud- 
ies. However, there are also IMF variations when observ- 
ing, analyzing, and interpreting the dat a of several clus- 
ters in the same manner. For example, iPhelps fc Janes! 
(11993') derive the IMFs of eight young open clusters in a 
consistent way. Although the average value of the power 
law index F = — 1.4±0.13 is in agreement with the 
Salpeter IMF, two of the clusters show strong deviations 

with r ~ - 1.8 (NGC 581) and F 1.1 (NGC 663). 

IScalol (12005") mentioned that these observed variations in 
cluster IMFs might still be explained by the uncertainties 
and, if that is the case, then the cluster data are consis- 
tent with a universal IMF but with sizeable variations, 
preventing us from determining an average IMF or F. 

How does our resulting IMF of NGC 3603 w i th F = 
-0.74 fit with the fundamental hmit in lKroupal (|200in ? 
iKroupal (|2001h explains the observed scatter is mainly by 
the following three uncertainties: (1) Poisson noise due 
to the finite number of stars, (2) dynamical evolution 
of the star cluster, and (3) unresolved binaries. In our 
determination the Poisson noise is taken into account 
when fitting the power law. It is only AF ~ ±0.02 for 
our ~ 10, 000 stars. As discussed in § [H we estimate 
that the dynamical evolution in NGC 3603 is expected to 
be insignificant for the intermediate- and low-mass stars 
and thus does not affect significantly the determination 
of the IMF from the field covering up to r < 110". This 
is supported by the fact that the IMF does not steepen 
significantly b eyond r > 30". As for unresolved binaries, 
iKroiipal (HooT) der ives that the correction for unresolved 
binaries will typically steepen the IMF by 0.5 < AF < 
0.8 for stars within 0.08 — 1 Mq, but it does not change 
significantly for > 1 Mq . Our simulation of the potential 
impact of unresolved binaries on the IMF within 0.4 — 20 
Mq yields AF ^ —0.04 (see ^ 19. 9p . Since the mass range 
is diffe rent in our analysis and the simulation of Kroupa 
(|2001h . it should not directly be compared. However, 
it is worth mentioning that there is a difference in the 
adopted binary frequencies. From the studies of the ONC 
we adopt a binary frequency of 15 and 30% for 0.1 — 2 and 
> 2 Mq, respectively, while IKroupal (|2001[ ) adopts 100% 
for > 3 Mq. Our result is similar to that of lStolte et al.l 
(|2006[ ). in which their correction of unresolved binaries 
steepens, by AF ~ -0.06, the IMF for stars within 0.4 - 
20 Mq in the field 7" < r < 33". 

We conclude that our analysis of the three main uncer- 
tainties find the IMF of NGC 3603 to be distinguishably 
flatter than the Salp eter IMF. But it is still within the 
fundamental limit in IKroupal ( 20011). Theref ore, as long 
as the intrinsic uncertainty by Kroupal (|200lD accommo- 
dates an IMF with a power law index of F = —0.74, 
we are not yet able to give a conclusive answer to the 
question of the IMF variation. 

10.3. Is the IMF top heavy in young and massive 
starburst clusters? 

The above IMF studies suggest that the IMF is al- 
most universal in many different star-forming environ- 
ments but still with nonnegligible variations in the power 
law index. Although there is no clear systematic trends 
in these variations, there is some evidence that these vari- 
ations are related to the stellar density. Compared to the 
field IMF, the IMF for intermediate- and high-mass stars 
tends to be slightly steeper in sparsely populated star- 



forming regions, while it is sometimes slightly flatter in 
very dense star clusters. 

Among young massive star-forming clusters, we find 
NGC 3603, the Arches cluster, and the GC cluster with 
a slightly fiatter IMF. For the Arches cluster iFiger et aTl 
fl999.) derive an IMF power la w index of F ^ —0.6 5 for 
stellar masses down to 10 Mq. IStolte et al.l (|2005f ) find 
the present-da y mass func t ion (P DMF) with F = —0.86 
for 6 - 60 Mq. iKim et al.l ()2OO60 derive the PDMF with 
F ^ —0.91 for 1.3 — 50 Mq. Using numerical simula- 
tions, they also correct for the dynamical evolution of 
the cluster and trace back the IMF to have a power 
law index F ^ —1.0 to —1.1. Interestingly, their PDMF 
has a flatter slope (F = —0.7 1) if only high-m ass stars 
are considered (5 — 50 Mq). iDib et ah (2007) recently 
reproduce this shallow IMF in the high-mass range in 
the Arches cluster in their models investigating the co- 
alescence and subsequent collapse of pre-stellar cores in 
molecular clouds. This suggests a primordial origin of the 
IMF in the Arches cluster. Also the young stellar popu- 
lation in the GC cluster exhibits a top heavy IMF. From 
the observed if-band LF of the high-mass stars orbiting 
the c entral massive black h ole in two counter-rotating 
disks. iPaumard et al.l(|2006f) report that the IMF is likely 
substantially flatter (by 1—1.5 dex) than the Salpeter 
IMF. However, these stars probably have formed in very 
dense gas disks rather than in a self-contracting cloud, 
and thus a direct comparison with NGC 3603 and other 
clusters is somewhat problematic. 

A common characteristic of these clusters, which al- 
ways show a slightly fiatter IMF, is the high central mass 
density. The core mass density of NGC 3603 is at least 
'--^ 6 X 10^ Mq pc~^. Note that this value is derived from 
the volume mass density profile, and the central density 
directly measured from the surface density (as sometimes 
given in literature for other clusters) is a factor of ~ 2 
larger (see 18.21) . The Arches cluster has a total mass of 



1 X 10* Mq similar to that of NGC 3603 and a core 
mass density 3 x lO'^ Mq pc"^ (|Figer et al.lll999[ ). sev- 
eral times larger than NGC 3603. The stellar density in 
the core of the GC cluster is expected to be much larger 
due to the presence of the central massive black hole. For 
example the stellar density in the central ar csecond of the 
GC i s estimated to be > 3 x 10'' Mq pc~^ (|Genzel et al.l 
I200I . 

In contrast, the central mass densities in the other 
Galactic or Local Group young stellar clusters are about 
an order of magnitude smaller than those in NGC 
3603 and the Arches cluster. The ONC, for exam- 
ple, has a core radius of ~ 0.2 pc, which is similar 
to NGC 3603, but has an order of magnitude smaller 
total mass of ^ 10^ Mq. It has thus an order of 
magnitude smaller central mass density of 2 x 10* Mq 
pc-3 (|Hillenbrand fc Hartma^ll998D . The IMF of the 
Trapezium cluster i n the ONC follows t he Salpeter IMF 
above a solar mass (jMuench et al.|[2"002[ ). 

An example of a cluster that does not exactly follow 
this line of argument is the R136 cluster in 30 Dor. Hav- 
ing a central mass density comparable to those in NGC 
3603 and the Arches cluster, the IMF of R136 is only 
marginally different from a Salpeter-like IMF. But there 
is still a considerable variation among different studies. 
Hosting mo re than 50 03- type, WN 6, and 0/WN6 tran- 
sition stars (|Massev fc Hunteilll998[ ) and, having a total 
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mass of 2 — 3 X 10^ Mq, R136 has been estimated to 
have a central mass density of ~ 5.5 x 10'' Mq pc^^ 
from stars > 2.8 Mq within r < 0.11 pc, which coul d 
be 3 times larger f or > 0.1 Mq ()Hunter et al.lll995D . 
iMackev fc Gilmord (|2003f ) de rive 3x 10^ M (t , pc"^ within 
r < 0.32 pc. As for its IMF, iBrandl et all ()1996f ) derive 
power law indices of F = —1.3 at r < 0.4 pc and F = —2.2 
at r > 0.8 pc, and find an average IMF wit h F = —1.6 
for th e high-mass stellar population. Massev fc Hunte"^ 
(|1998f) measure the IMF with F = -1.3 to - 1.4 for 
2.8 - 120 Mq. In contrast. [Sirianni et all ()2000[) derive 
a broken power law IMF with F = —1.28 for 2.1 — 6.5 
Mq and F = -0.27 for 1.35 - 2.1 Mq, indicating a sub- 
stantially flat IMF over the entire mass range. Although 
there is ample evidence for a flattening of the IMF at 
subsolar masses, the observed flattening at such high 
masses of ^ 2 Mq in this study is somewhat peculiar. 
Thus it might be the case that R136 is in conflict with 
our hypothesis of the flat trend of the IMF in the dense 
star clusters. It is not yet possible for us to argue in 
more detail about the IMF of R136, as a reliable mea- 
surement for low- and intermediate-mass stars is more 
difficult because of its much larger distance than those 
of the Galactic clusters. 

As a summary, although with low statistical signifi- 
cance, the variations of the IMF could be linked to the 
stellar density of the population. In particular the most 
massive starburst clusters like NGC 3603 and the Arches 
cluster with their dense central regions tend to have a 
distinguishably flatter IMF compared to the standard 
Salpeter-hke IMF. 

11. SUMMARY AND CONCLUDING REMARKS 

Our study is aiming at a fundamental question of cur- 
rent star formation research. Is the IMF universal, or 
does it vary with environment? To answer this question, 
we measured the IMF of NGC 3603 - one of the most 
massive Galactic star-forming regions - from our NIR ob- 
servations with the AO system NACO at the VLT/ESO. 
In the following we summarize the main results from our 
study. 

1. NIR Photometry 

Our very deep, high angular resolution JHKgL'- 
band images obtained by NACO show unprece- 
dented details of the core of the starburst cluster in 
NGC 3603. Together with the wider field ISAAC 
JHKs images, we could successfully derive mag- 
nitudes and positions of almost 10,000 stars in the 
dense cluster up to r ~ 110" covering the mass 
range from the most massive stars down to ~ 0.4 
Mq. The brightest 256 stars in the NACO images 
of the inner r < 13" are listed and cross-identified 
with potential counterparts in previous studies. 

2. Age, extinction, and disk fraction 

Based on the fitting of the stellar evolution models 
in the CMDs and CCDs, we derived the age of the 
PMS stars of 0.5 - 1.0 Myr and the upper limit 
for the age of the MS stars of ~ 2.5 Myr, suggest- 
ing a slight age spread in the cluster. The derived 
average foreground extinction is Ay = 4.5 ± 0.5 
mag, and the foreground extinction increases by 
AAy 2.0 mag towards larger radii (r > 55"). 



Using the Ks — L' versus J ~ H CCD, we derived a 
circumstellar disk fraction of ~ 25 ± 10% for stars 
with a mass of > 0.9 Mq in the central cluster 
(r < 10"). 

3. IMF and its radial variation 

Applying the field star rejection and the incom- 
pleteness correction, the KIjF for 7514 stars si- 
multaneously detected in the JHKs bands follows 
a power law with no obvious turnover or trunca- 
tion within the detection limit oirriKs ^ 17.4 mag 
(based on the J-band 50% completeness of ~ 19.4 
and the typical color of J — Ks ^ 2). Within the 
mass range of 0.4 - 20 M© the IMF is well de- 
scribed by a single power law with a power law 
index F ~ —0.74. We found the power law index 
decreasing from F ~ —0.3 at r < 5" to F ~ —0.8 at 
r ~ 30". The strong steepening occurs in the inner 
r < 13", pointing towards mass segregation in the 
very center of the cluster. No significant variation 
of the IMF is found for larger radii (r > 30"). 

4. Size, mass, and dynamical status 

Fitting a King model to the radial density profile 
of stars with a mass of 0.5 — 2.5 Mq, we derived 
a core radius of ~ 4". 8 (~ 0.14 pc at d ~ 6 kpc). 
As the radial density decreases even at the limits 
of our field of view, we can give a firm lower limit 
of r = 110" (~ 3.2 pc) for the cluster size. We also 
derive an upper hmit of r = 1260" 37 pc) for 
the tidal radius of the cluster. The de-projected 
King model allowed us to extrapolate to the total 
mass of NGC 3603. Assuming a single power law 
IMF with index F = —0.74 within the mass range 
of 0.1 — 100 Mq, we found a total mass of about 
1.0 - 1.6 X 10-* Mq. The half-mass radius is found 
to be within 25" - 50" (0.7 - 1.5 pc). The derived 
core mass density of the cluster is > 6 x 10** Mq 
pc~^. 

We estimate a half-mass relaxation time of approx- 
imately 10 — 40 Myr for stars with a typical mass 
of 1 Mq , an order of magnitude larger than the age 
of the PMS population in the cluster (< 1 Myr). 
This implies that the intermediate- and low-mass 
stars have not yet experienced significant dynam- 
ical relaxation. However, the relaxation time of 
the high-mass stars is expected to be an order of 
magnitude shorter and is comparable to the clus- 
ter age. We could thus not conclude if the observed 
mass segregation of the high-mass stars is caused 
by dynamical evolution or if it is primordial. In- 
deed it can be due to the combination of them. We 
compute that our images with a maximum radius 
of r ~ 110" cover at least ~ 67% of intermediate- 
and low-mass stars of NGC 3603. The stars out- 
side the observed field cannot steepen the IMF by 
more than AF < 0.16. Considering also the fact 
that the IMF does not significantly change beyond 
r > 30", we conclude that the observed IMF is rep- 
resentative for the whole NGC 3603 stellar cluster, 
irrespective of the mass segregation in the very cen- 
ter. 

5. Systematic uncertainties of the IMF 

We thoroughly analyzed the systematic errors in 
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the IMF determination. In particular we derived 
the errors from uncertainties in the age, distance, 
foreground extinction, stellar evolution models, 
metallicity, individual extinction, incompleteness 
correction, cluster membership, unresolved bina- 
ries, and the stellar mass outside the observed field. 
Combining all errors we derived the power law in- 
dex of r = — 0.741q;47. Assuming a Gaussian prob- 
ability distribution, we conclude that the probabil- 
ity that the IMF is as steep as the Salpeter IMF 
(r = -1.35) is less than - 10%. 

Our result thus supports the hypothesis of a top- 
heavy IMF in massive star-forming clusters and starburst 
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galaxies. One potentially common property among such 
clusters showing flat IMFs would be the high stellar den- 
sity in the core of the cluster. 
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APPENDIX 
PHOTOMETRY LIST 

Thanks to our AO-assisted NACO data with the unprecedented high resolution, we successfully identify objects in 
the highly crowded cluster center and perform photometry from the brightest WN stars to the faintest stars in subsolar 
masses. Here we present a photometry list (Table 6) of JHKs{L')-hand detected stars for magnitude mj < 15.5 mag 
(> 95% completeness) in the inner r < 13" region (diameter of ~ 0.75 pc at d ~ 6 kpc). This list encompasses all 
sources within the magnitude range on the projected plane, and any rejections of field stars have not been performed. 

The list contains 256 detected sources. Each source has an ID number, a coordinate, JHKs magnitudes, and a 
L' magnitude in the case of a source detected in the all four bands. The sources are aligned according to increasing 
mj. The coordinates are indicated in their offsets from the cluster center [RA = ll''15'"05^90, Dec = -61°15'28".4 
(J2000)]. We note that, since the photometric calibration of the L'-band photometry was done by merely one data set 
of a standard star, the presente d L'-band magnitude has to be considered as a roughly calibrated photometry. The 
spectral types are adopted from ICrowther fc DessartI (Il998l Tab l e 8 th erein), which consists of their stellar spectral 
types combined with those from lMoffatI (|1983l ) and lDrissen et al.l ()1995D . 

In case of sources which are identical or likely to be counterparts to objects in other studies, the design ations or IDs 
are pr esented. We perform cross-identifications of sources with the HST optical detections presented in lMoffat et al.l 
(|l994f ) (MDS94), with the ground-based mid-IR detections in ^iir nbcrgcr & Stankc (2003) (NS03), and with the 
optical detections {V < 14 mag) from the HST data set in ISung fc BesselH ((2004) (SB04). The identifications of 
sources are done by means of a visual inspection after making maps of the examining lists. We obtain the impression 
that there seems to be some degree of position mismatches across fields among the maps due perhaps to differences of 
instrumental properties such as resolutions, sizes of field of view, and geometric distortions. Thus, we chose to use the 
eye inspection rather than applying an y solid criteria , e.g., a tolerant distance of a coincidence. As for designations 
of the counterparts, the source IDs in iMelnick et al.l (I1989D fMTT), a s well as the letters Al-3 to F (jvan den BosI 
Il928t iHofmann fc WeigelH 1 1986( 1 are simply adopted from iMoffat et al.l (|1994[ ). In addition, the X-ray counterparts 
obtained by the Chandra satelli te in iMoffat et al.l (120021) are indicated. We simply adopt the MDS94 sources with 
X-ray counterparts presented in IMoffat et al.l (|200^ . Figure 2 therein) . Only bright WR and early 0-type stars with 
X-ray counterparts are presented, and since the authors suggested that many detected X-ray sources might be low- 
mass PMS members, many more sources in our NACO list are expected to ha ve X-ray counterparts . As t he spatial 
scale of our NACO NIR data is much smaller than that of the mid-IR data in iNiirnberger fc Stankd ()2003l ). a single 
mid-IR source could be superposed with several NIR sources. In such a case, the source that would have a dominant 
flux contribution in the NIR is presented. 



24 



Harayama et al. 



REFERENCES 



Armitage P. J., Clarke C. J., Palla F. 2003, MNRAS, 342, 1139 
Ascenso, J., Alves, J., Beletsky, Y., Lago, M. T. V. T. 2007, A&A, 
466, 137 

Balick, B., Boeshaar, G. O., & Gull, T. R. 1980, ApJ, 242, 584 
Baraffe, I., Chabrier, G., AUard, F., & Hauschildt, P. H. 1998, 
A&A, 337, 403 

Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 2002, 

A&A, 382, 563 , 

Barlow, R. J. 2003, preprint farXiv:physics/03_06168 ) 

Bertelli, G., Bressan, A., Chiosi, 0., f agotto, t ., & Nasi, E. 1994, 

A&AS, 106, 275 
Bessell, M. S., & Brett, J. M. 1988, PASP, 100, 1134 
Binney, J., Tremaine, S. 1987, Galactic Dynamics, Princeton 

University Press, 514 
Bonatto, Ch., Bica, E., & Girardi, L. 2004, A&A, 415, 571 
Bonnell, I. A., & Davies, M. B. 1998, MNRAS, 295, 691 
Bouvier, J., Stauffer, J. R., Martin, E. L., Barrado y Navascues, 

D. , Wallace, B., & Bejar, V. J. S. 1998, A&A, 336, 490 
Brandl, B., Sams, B. J., Bertoldi, F., et al. 1996, ApJ, 466, 254 
Brandl, B., Brandner, W., Eisenhauer, F., Moffat, A. F. J., Palla, 

F., & Zinnecker, H. 1999, A&A, 352, L69 
Brandner, W., Grebel, E. K., Chu, Y.-H., & Weis, K. 1997, ApJ, 
475, L45 

Bricefio, C, Luhman, K. L., Hartmann, L., Stauffer, J. R., & 

Kirkpatrick, D. 2002, ApJ, 580, 317 
Brown, A. G. A. 1998, ASP Conference Series, 142, ed. G. Gilmore 

& D. Howell (San Francisco), 45 
Clayton, C. A. 1986, MNRAS, 219, 895 
Crowther, P. A., & Dessart, L. 1998, MNRAS, 296, 622 
Crowther, P. A., Lennon. D. J.. Walborn. N. R., & Smartt, S. J. 

2006, preprint ( ar_Xiv;astro-ph/0606717 ) 
Chabrier, G. 2003, rAsr, ilo, I06 

Demarque, P., Woo, J.-H., Kim, Y.-C, & Yi, S. K. 2004, ApJS, 
155, 667 

De Free, C. G., Nysewander, M. C, & Goss, W. M. 1999, AJ, 117, 
2902 

Dib, S., Kim, J., & Shadmehri, M. 2007, MNRAS, 381, L40 
Diolaiti, E., IJendinelli, O., Bonaccini, D., Close, L., Currie, D., & 

Parmeggiani, G. 2000, A&A, 147, 335 
Djorgovski, S., et al. 1995, ApJ, 438, L13 

Drissen, L., Moffat, A. F. J., Walborn, N. R., & Shara, M. M. 1995, 
AJ 110 2235 

Duchlne, G., Bouvier, J., & Simon, T. 1999, A&A, 343, 831 
Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485 
Eisenhauer, F., Quirrenbach, A., Zirmecker, H., & Genzel, R. 1998, 
ApJ, 498, 278 

Elson R. A. W., Fall M. S., & Freeman K. C. 1987, ApJ, 323, 54 
Figer, D. F., Kim, S. S., Morris, M., Serabyn, E., Rich, R. M., & 

McLean, I. S. 1999, ApJ, 525, 750 
Figueredo, E., Blum, R. D., DamineU, A., & Conti, P. S. 2002, AJ, 

124 2739 

Fischer, D. A., & Marcy, G. W. 1992, ApJ, 396, 178 

Frogel, J. A., Persson, S. E., & Aaronson, M. 1977, ApJ, 213, 723 

Genzel, R., et al. 2003, ApJ, 594, 812 

Girardi, L., Bressan, A., BerteUi, G., & Chiosi, C. 2000, A&AS, 
141, 371 

Goss, W. M., & Radhakrishnan, V. 1969, Astrophys. Lett., 4, 199 
Grabelsky, D. A., Cohen, R. S., Bronfman, L., & Thaddcus, P. 

1988, ApJ, 331, 181 
Grebel, E. K. 2004, ASP Conference Series, 322, ed. Lamers et al. 

(San Francisco), 101 
Hernandez, J., et al. 2007, ApJ, 662, 1067 

Hillenbrand, L. A., Strom, S. E., Calvet, N., Merrill, K. M., Gatley, 
I., Makidon, R. B., Meyer, M. R., & Skrutskie, M. F. 1998, AJ, 
116, 1816 

Hillenbrand, L. A., & Hartmann, L. W. 1998, ApJ, 492, 540 
Hillenbrand, L. A., & White, R. J. 2004, ApJ, 604, 741 
Hofmann, K.-H., Weigelt, G. 1986, A&A, 167, L15 
Hofmann, K.-H., Seggewiss, W., Weigelt, G. 1995, A&A, 300, 403 
Hunter, D. A., Shaya, E. J., Holtzman, J. A., Light, R. M., O'Neil, 

E. J., Jr., & Lynds, R. 1995, ApJ, 448, 179 
Kennicutt, R. C, Jr. 1984, ApJ, 287, 116 
Kenyon, S. J., & Hartmann, L. 1995, ApJS, 101, 117 
Kenyon, S. J., & Gomez, M. 2001, AJ, 121, 2673 
Kerr, F. J. & Lynden-Bell, D. 1986, MNRAS, 221, 1023 

Kim, S. S., Figer, D. F., Kudritzki, R. P., & Najarro, F. 2006, ApJ, 

653, L113 
King, I. 1962, AJ, 67, 471 
Kroupa, P. 2001, MNRAS, 322, 231 
Kroupa, P. 2002, Science, 295, 82 
Kohler, R. 2001, AJ, 122, 3325 

Kohler, R., Petr-Gotzens, M. G., McCaughrean, M. J., et al. 2006, 
A&A, 458, 461 

Kouwenhoven, M. B. N., Brown, A. G. A., Zinnecker, H., Kaper, 

L., & Portegies Zwart, S. F. 2005, A&A, 430, 137 
Lada, C. J., & Adams, F. 1992, ApJ, 393, 278 



Lada, C. J., Muench, A. A., Haisch, K. E., et al. 2000, AJ, 120. 
3162 

Lejeune, T., & Schaerer, D. 2001, A&A, 366, 538 

Luhman, K. L., Stauffer, J. R., Muench, A. A., Rieke, G. H., Lada, 

E. A., Bouvier, J., & Lada, C. J. 2003, ApJ, 593, 1093 
Mackey, A. D., & Gilmore, G. F. 2003, MNRAS, 338, 85 
Maercker, M., & Burton, M. G. 2005, A&A, 438, 663 
Maercker, M., Burton, M. G., & Wright, C, M. 2006, A&A, 450, 

253 

Malkov, O., Zinnecker, H. 2001, MNRAS, 321, 149 

Massey, P., & Hunter, D. A. 1998, ApJ, 493, 180 

Melnick, J., & Grosb0l, P. 1982, A&A, 107, 23 

Melnick, J., Tapia, M., & Terlevich, R. 1989, A&A, 213, 89 

Meynet, G., & Maeder, A. 2003, A&A, 404, 975 

Miller G. E., Scalo J. M. 1979, ApJS, 41, 513 

Moffat, A. F. J. 1983, A&A, 124, 273 

Moffat, A. F. J., Drissen, L., & Shara, M. M. 1994, ApJ, 436, 183 
Moffat, A. F. J., et al. 2002, ApJ, 573, 191 

Moffat, A. F. J., Poitras, V., Marchenko, S. V., Shara, M. M., 
Zurek, D. R., Bergeron, E., Antokhina, E. A. 2004, AJ, 128, 
2854 

Muench, A. A., Lada, E. A., Lada, C. J., Alves J. 2002, ApJ, 573, 
366 

Niirnberger, D. E. A., & Petr-Gotzens, M. G. 2002, A&A, 382, 537 
Niirnberger, D. E. A., Bronfman, L., Yorke, H. W., & Zinnecker, 

H. 2002, A&A, 394, 253 
Niirnberger, D. E. A., & Stanke, T. 2003, A&A, 400, 223 
Palla, F., & Stahler, S. W. 1999, ApJ, 525, 772 
Pandey, A. K., Ogura, K., & Sekiguchi, K. 2000, PASJ, 52, 847 
Patience, J., Chez, A. M., Reid, I. N., Weinberger, A. J., Matthews, 

K. 1998, AJ, 115, 1972 
Paumard, T., et al. 2006, ApJ, 643, 1011 

Petr-Gotzens, M. G., Coude du Foresto, V., Beckwith, S. V. W., 

Richichi, A., & McCaughrean, M. J. 1998, ApJ, 500, 825 
Pfalzner, S., Olczak, C, Eckart, A. 2006, A&A, 454, 811 
Phelps, R. L., & Janes, K. A. 1993, AJ, 106, 1870 
Pinfield, D. J., Hodgkin, S. T., & Jameson, R. F. 1998, MNRAS, 
299, 955 

Portegies Zwart, S. F., Makino, J., McMillan, S. L. W., & Hut, P. 

2002, ApJ, 565, 265 
Preibisch, Th., Balega, Y., Hofmann, K.-H., Weigelt, G., Zinnecker, 

H. 1999, New Astronomy, 4, 531 
Preibisch, T., Stanke, T., & Zinnecker, H. 2003, A&A, 409, 147 
Prisinzano, L., Micela, G., Sciortino, S., & Favata, F. 2003, A&A, 

404, 927 

Prosser, C. F., Stauffer, J. R., Hartmann, L., et al. 1994, ApJ, 421, 
517 

Reid, L N., & Gizis, J. E. 1997, AJ, 113, 2246 

Rieke, G. H., & Lebofsky, M. J. 1985, ApJ, 288, 618 

Rudolph, A. L., Fich, M., Bell, G. R., Norsen, T., Simpson, J., P., 

Haas, M. R. & Erickson, E. F 2006, ApJS, 162, 346 
Sagar, R., Munari, U., & de Boer, K. S. 2001, MNRAS, 327,23 
Sapleter, E. E. 1955, ApJ, 121, 161 

Sanner, J., Altmann, M., Brunzendorf, J., & Geffert, M. 2000. 

A&A, 357, 471 
Scalo, J. M. 1986, Fundamentals of Cosmic Physics, 11, 1 
Scalo, J. M. 1998, in ASP Conf. Ser. 142, The Stellar Initial Mass 

Function (38yh Herstmonceux Conference), cd. G. Gilmore & D. 

Howell (San Francisco: Astron. Soc. Pac), 201 
Scalo, J. M. 2005, in IMF@50: The Initial Mass Function 50 Years 

Later, ed. E. Corbelli, F. Palla, & H. Zinnecker (Dordrecht: 

Kluwer), 23 

Schmutz, W., Drissen, L. 1999, Rev. Mexicana Astron. Astrofis. 

Ser. Conf., 8, 41 
Shatsky, N., & Tokovinin, A. 2002, A&A, 382, 92 
Sher, D. 1965, MNRAS, 129, 237 

Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593 
Sirianni, M., Nota, A., De Marchi, G., Leitherer, C, & Clampin, 

M. 2000, ApJ, 533, 203 
Slesnick, C. L., Hillenbrand, L. A., & Massey, P. 2002, ApJ, 576, 

880 

Sollima. A.. Ferraro. F. R., & Bellazzini, M. 2007, 

arXiv: astro-ph /0708072 
Stolte, A., arandner, W., Brandl, B., Zinnecker, H., & Grebel, E. 

K. 2004, AJ, 128, 765 
Stolte, A., Brandner, W., Grebel, E. K., Lenzen, R., & Lagrange, 

A.-M. 2005, ApJ, 628, L113 
Stolte, A., Brandner, W., Brandl, B., Zinnecker, H., & Grebel, E. 

K. 2006, AJ, 132, 253 
Sung, H., & Bessell, M. S. 2004, AJ, 127, 1014 
Tapia, M., Bohigas, J., Perez, B., Roth, M., Ruiz, Ma. T. 2001, 

RMxAA, 37, 39 
Tout, C, Livio, M., & Bonnell, I. 1999, MNRAS, 310, 360 
van den Bergh, S. 1978, A&A, 63, 275 

van den Bos, W. H. 1928, Bull. Astron. Inst. Netherlands, 4, 261 
Walborn, N. R. 1973, ApJ, 182, L21 



The IMF of NGC 3603 



25 



Walborn, N. R., Howarth, I. D., Lcnnon, D. J., et al. 2002, AJ, 
123, 2754 

Wilking, B. A., Meyer, M. R., Greene, T. P., Mikhail, A., & 

Carlson, G. 2004, AJ, 127, 1131 
Yan, L., McCarthy, P. J., Storrie-Lombardi, L. J., & Weymann, R. 

J. 1998, ApJ, 503, L19 



Yi, S., Demarque, P., Kim, Y. -C, Lee, Y.-W., Ree, C. H., Lejeune, 
T., & Barnes, S. 2001, ApJS, 136, 417 



TABLE 6 Photometry list 



Num. ARXra ADecTW" 



H 



Ks L' Type'fai 



"NSO3P] SBO4H 



0.53 



-7.15 



10.28 



7.78 
7.98 
8.49 
9.38 

9.68 
9.94 
9.95 
10.09 
10.16 
10.24 
10.25 
10.29 



m 



10.71 



11.67 10.82 10.71 

11.68 11.28 11.16 
11.76 11.36 11.13 
11.79 11.27 11.25 
11.81 11.36 11.19 
11.83 12.90 13.04 
11.88 11.59 11.44 

.2.05 12.62 



11.43 
11.67 



WN6h-|-abs 
WN6h-|-abs 
WN6h+abs 

03 V 

03 III 
05 Ill(f) 

03 III , ^ 
05.5 Ill(f) 

05 V 

04 V 

Q3 V, , 

04 V(f) 

05 V 

il 

04 V 

04 V 

05 V 




^-III 
M V 

36.5 V+? 
04 V 
04 V 



04 V 
04 V 



12.10 11.56 



1.99 



H23(B) 

Hl<=9l30(Al) 

H18(C) 

H['^9l31(A2) 

H''^''l26{A3) 

H22(17) 

H42^ , ^ 

H39 6/F) 

H19(E) 

Hlciol49(DbB]) 

H40, , 
H51(23) 

Hl<=iol5o(D['^sl) 
H16 

25 

38 , , 
H62(15['^^l) 
H9{2lM) 
H6l''=il(10/G) 

S'm , , 
58(15) ["=31 

H7(26) 



43 , , 

52(D[<:5]) 



(H7b) 10(22 [='^1) 
68(44['^61) 

8^21 M) 
12(30''^81) 

H6l['^il(10/G) 
34 

70 , , 
13(30l'=8]) 

ll(22['^''l) 



6N[d31 
6K[d21 

6G[dil 
6N[d31 



6G[diI 



6K[d21 



6M 



10582 
10558 
10635 
10559 
10570 

10478 
10392 
10543 
10422 
10498 
10616 

10423 
10706 
10534 
10419 
10396 
10895 
10207'°y 

10683 



10363 
10705 



10319 
10361 



10207[«y 



11.66 
11.58 
11.82 
11.85 



12.28 



11.50 
11.63 
11.57 
11.29 

IIM 

11.61 
11.69 



.42 



.71 



21, , 
67[<=2](36) 

67['^2](36) 

53(D[<=5]) 

54 , , 
69(44l<=6]) 



6G[dil 



04 V 



29 



60 



26 
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Num. AR.A, 



ADecl 



H 



Ks 



L' Type 



MDS94l'^ 



Note. — Photometry and counterparts of 256 JHKs{L )-deteeted sourees brighter than mj — 15.5 (~95% eompletencss) in the central region 
r < 13" of the star-forming cluster in NGC 3003. 

^'^^ Source location is shown in an offset from the center of the cluster (RA — Ilhl5m05.90s, Dec — — 61°15'28.4" [J2000]). 

Spectral type is a dopted from a lis t in [Cro^^er i^^Des sartl (119981 . Table 8 therein), which consists of spectral types in their study combined 
with those in [Moffat 1 1983) and in lDrissen et al.l l ll995ir 

Designated ID s of potential counte rparts in HST optical detections bv lMoff at et ar|)|1994| ) . A so urce with prefix H is found to have a X-ray 
counterpart bv'Moffat et al.l ll2002D . IDs originated in previous wo rks such as[M e^^^ ^et^ll jl989f) (in parentheses), and letters were adopted 
from tables in Moffat et al. il994D. The letters A-F originated in Ivan den Boj I ll928i) (G in Walborn 1973), and the source A was resolved 
into Al-3 bv .Hofmann fc Weigeld II1986I '). 

'■=^1 Star 21 and 52 together would correspond to MDS-61 (G/MTT-10). 

'■=^1 Star 73 and 74 together would correspond to MDS-67 (MTT-36). 

'■=^1 Multiple detections. MDS-58 and 62 correspond to MTT-15. 

MDS-8 and 9 correspond to MTT-21. 
'■^^l MDS-49, 50, 52, and 53 correspond to D. 

MDS-68 and 69 correspond to MTT-44. 

MDS-10 and 11 correspond to MTT-22. 

MDS-12 and 13 correspond to MTT-30. 
[c9] j^-^^ and A3 together consist of a X-ray source within the match of the smoothed PSF in lMofTat et all j2002D . 
[clo] ]yi]3g_49 ^j^j 50 together consist of a X-ray source. 



Designated IDs of potential counterparts in mid-infrared detections by INiirnberger Sz Stankd 1 120031) . The 11.9 /im sources 6A-F are in the 
concentrated cluster core. 

f^^^ Star 19, 33 and 51 together would correspond to 6G. 

^^^^ Star 18 and 43 together would correspond to 6K. 

I'*^' Star 17 and 22 together would correspond to 6N. 



IDs in lSung fc Besselll I I2004I . Table 2 therein). Sources brighter than V — 14 mag detected in their analysis of the HST PCI archival data 
set. 



f^"*"^ Star 21 and 52 together would correspond to 10207. 



